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ABSTRACT 

While star-forming galaxies could be major contributors to the cosmic GeV 7-ray background, they are 
expected to be MeV-dim because of the "pion bump" falling off below ^100 MeV. However, there are very few 
observations of galaxies in the MeV range and other emission processes could be present. We investigate the 
MeV background from star-forming galaxies by running one-zone models of cosmic ray populations, taking 
into account the leptonic emission, including Inverse Compton and bremsstrahlung, as well as nuclear lines 
(including ^''Al), emission from core collapse supemovae, and positron annihilation emission, in addition to 
the pionic emission. We use the Milky Way and the GeV-TeV detected starbursts M82 and NGC 253 as 
templates of normal and starburst galaxies, and compare our models to radio and GeV-TeV 7-ray data. We 
find that (1) Inverse Compton losses off the CMB flatten out the pion bump at high z for normal galaxies, 
(2) we cannot rule out that starbursts have significant MeV emission if their magnetic field strengths are low, 
and (3) cascades can contribute to the MeV emission of starbursts if they emit mainly hadronic 7-rays. The 
star-forming galaxy contribution to the GeV background is uncertain by an order of magnitude, depending 
on how much of the cosmic star-formation is in starbursts. Our fiducial model predicts that ~ 1/3 of the 
unresolved GeV background is from star-forming galaxies, with comparable contributions from normal and 
starburst galaxies. About ^ 2% of the claimed 1 MeV background is diffuse emission from star-forming 
galaxies; we place a firm upper limit of ^ 10% contribution based on the requirement that star-forming galaxies 
do not overpower the observed 7-ray background at any energy. The low star-forming galaxy contribution 
arises because the observed 7-ray background spectrum is steeply falling with energy, while the star-forming 
contribution is slowly increasing with energy in the MeV range. A different source class must emit the observed 
MeV background, if it is real. 

Subject headings: gamma rays: diffuse background - gamma rays: galaxies - cosmic rays - galaxies: starburst 



1. INTRODUCTION 

The origin of the cosmic 7-ray background (Figure E} re- 
mains a mystery. The unresolved extragalactic 7-ray back- 
ground appears to have a single power law spectrum ex- 
tending from ^ 50 Me V to 100 GeV, as determined from 
EGRET and Fermi-LAT ( SreekumaretaDl 19981: IStrong et all 
|2()04a; Abdo et al. 2010a). There has been much recent de- 
bate about the origin of this emission, particularly whether 
it comes from blazars or star-forming galaxies. The calcu- 
lated blazar contribution (dominated by the soft Flat Spec- 
trum Radio Quasars at low energies and the hard BL Lacs at 
high energies; Pavlidou & Venters 2008i IVenters & Pavlidoul 
I2OI ll) has a similar spectrum to the 7-ray background and, if 
they make up most of the unresolved backgro und, the remain- 
ing blazars should be resolved fairly easily (lAbazaiian et alj 
I2OI II) . Star-forming galaxies are expected to have a "bump" 
in their emission at around a GeV, because most of their 7-ray 
emission comes from pion production b y cosmic ray protons 
with a threshold of a few hundred MeV (IProdanovic & Fieldi 
l2004l) . In contrast to blazars, ver y few star-forming galax - 
ies will be resolved with Fermi jPavlidou & Fields! 1200 lb . 
The contribution of both source classes at GeV energies is 
still unclear. The MeV background potentially is a pow- 
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erful way of discriminating the two sources: the SED of 
blazars is expected to be higher at MeV energies, while the 
SED of star -forming galaxies is expected to be lower (e.g., 
IStecker & V enters 20 11). 

At GeV energies, Fermi has found that the unre- 
solved extragalac tic 7-ray backgro und is a power law 
dN/dE cx E-^ '^^ (lAbdo etal.1 l2010ah . Before Fermi, the 
GeV 7-ray background was usually attributed to blazars, 
which dominate the population of detected 7-ray sources 
(Padovani et al.'1993; Stecker et al.'1993';'Salamon & SteckeiJ 
| 1994; Stecker & Salamon 1996; Mucke & Pohl 2000iK 
I Abdo et al.l ( 1201 Obi) claim that the Fermi source counts are 
inconsistent with a blazar origin of the unresolved GeV 
background, although th is conclusion has been disputed 
dStecker & Venters 2011; compare with Singal et allboill 
and Malyshev & Hogg 2011,) . and some estimates find that 
blazars are a minority of the unresol ved GeV background 
(iDermed 120071; llnoue & Totaiiil l2009l) . Studies of blazars 
detected in the flux-limited Australia Telescope 20 GHz 
survey also find that these blazars contribute a minority of 
the GeV b ackground, based on an observed 7-ray-radio 
correlation (iGhirlanda et al.l 1201 1 1 and Ferrni imag e stack- 
ing of GeV undetected blazars dZhou et all 1201 Ih . Other 
proposed alternatives to a blazar origin for the GeV 7-ray 
backgroun d include dark r natter annihilation and decay 
(e.g., Bergstrom et"aLl 120011: JUllio et al.|l2002t iBertone et alJ 
20071) intergalactic shocks Cpeb & Waxman 2000; Miniat^ 
2002tlKeshet et a l."2003'). non-blazar AGNs (Active Galactic 
Nuclei; e.g.. | In oue 2011a; Teng et al. 2011), galaxy clusters 
(iDar & Shavivlll995i: iCola fr ancesco & Blasi 1998 ), cascades 
from ultrahigh energy cosmic rays (e.g.. Wdowczvk et al.l 
[T97I ICoppi & AharonianlfT997l: iMuraseet al.llSoTl) . cosmic 
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rays accelerated by type la supernovae (iLien & Fieldsll2012b . 
and Galactic foreground emissi on from millisecond pulsars 
jFauch er-Giguere & Loebll2010h . cosmic rays in the Galactic 
corona ( Feldmann et aL | 2012 ), and unsubtracted diffuse 
emission (Keshet et al. 2004; Maly shev & Ho gg 201 1). 

Star-forming galaxies are also expected to contribute to the 
GeV 7-ray background. Supernovae and other processes ac- 
celerate cosmic rays (CRs), particularly protons but also other 
nuclei and electrons; these CRs can then interact with the 
ambient gas, radiation, and magnetic fields to emit broad- 
band radiation. EGRET only detected the Milky Way and 
the LMC (Abdo et al. 2010c), but Fermi has also de tected 
the SMC ( Abdo et aL,2010d) . M31 (|Abdo et al."2010e), and 
the starburst galaxies M82 and NGC 253 (Abdo et al. 20101), 
as well as the starbu rst-Seyferts NGC 4945 and NGC 1068 
dAbdo e t al. 2010g; Lenain et al.ll2010l) . The TeV telescopes 
VERITAS (Very Energetic Radiation Imaging Telescope Ar- 
ray System) and HESS (High Energy S tereoscopic System) 
have also detected M82 and NGC 253 (lAcciari et al.ll2009[ 
lAceroetal.ll2009h . The GeV emission is thought to come 
from the inelastic collisions of cosmic ray protons with pro- 
tons in the interstellar medium (ISM), creating pions which in 
turn decay into 7-rays, neutrinos, and secondary electrons and 
positrons. Leptonic processes, particularly Inverse Compton 
(IC) and bremsstrahlung, may also contribute at lower ener- 
gies. The sheer number of normal star-forming galaxies in- 
dicates that they contribute to the 7-ray background, even 
though they are individually faint, with pTe-Fenni estimates 
placing the contribution anywhere from a few percent to most 
of the GeV background (S trong et al.lll976l;lLichti et al.ll 19781; 
iPavlidou & Fields 2002). 

A special subset of star-forming galaxies are the less nu- 
merous but individually luminous starburst galaxies, which 
have high star-formation rates and thus high CR energy 
densities. The high gas densities can enhance the pio- 
nic 7-ray emission from starbursts (Torres et al. 2004) so 
that they are a major component of the GeV 7 -rav back- 
ground (iThompson. Ouataert. & WaxmanI 120071: see also 
iLoeb & Waxmanll2006h . It is even possible that starbursts are 
"proton calorimeters", converting the power injected as CR 
protons into pionic losses (Pohl 1994; Loeb & Waxman 200S 
IThompson. Ouataert. & WaxmanI 120071: ILackietal.1 1201 Oh 
but the extent to which t his is true is disputed (lAcero et al.l 
I2OO9I: iLacki et al.l IMl lOhm & HintonI lIoTl) . a s is the 
fraction of star-form ation in calorimetric galaxies (' Steckerl 
12007; Tho mpson et a l. 2006b; Stecker & Venters 201 li Re- 
cent estimates of the star-formation contribution to the GeV 
7-ray backgrou nd find values ra nging fro m 10% to 50% 
(e.g.. Bhattach arva & SreekumaH 12 009: Fi elds et al.l 1201 Ol 
iMakiya et aLi201 ll; IStecker & Venters»201 ID . with starbursts 
alone making up 1% to 50% of t h e GeV 7-ray background 
dMakiva et alj 1201 U iLacki et alj 1201 ll; IStecker & VentersI 
HoHF 

Although Fermi continues to advance our understanding 
of the GeV 7-ray sky, very little is known about the MeV 
7-ray background. The available data on the MeV back- 
ground comes from S olar Maximum Mission (SMM; 0.3 - 
7 MeV; I Watanabe et al. 1999,). Imaging Compton Tele scope 
(COMPTEL; 1 - 20 MeV; I Weidenspointner et al.ll2000l) . and 
Energetic Gamma Ray Experiment Telescope (EGRET; 3 
- 20000 MeV; Sree kumar et aH [T998t IStrong et al.l l2004ah . 
Blazars may constitute most of the MeV background 
fBloemen et al. 19 95; ZdziarskillT996t iGiommi et al.ll2006l: 
lAjello et al.ll2009l) . Additionally, 7-ray line emission and 



the accompanying Compton-downscattered continuum emis- 
sion from type la supernovae is expected to contribute 
("Clavton & Silk"1969^.'cia vton & Wardll975l;lThe et al.fT99l 
Ruiz-Lapuente et al. 2001), but more recent estimates show 
that supernovae account for only ^ 10% of the MeV 
backgroun d (Iwabuchi & Kumagaili200lUSt rigari et al.ll2005l: 
lAhn etal.1 [2005; Horiuchi & Beacom 2010). Other pro- 
posed sources include nonthermal electrons in AGN coronae 
(Inoue et al. 2008), the cores of misaligned ra dio loud AGNs 
(Inoue 2011a), the radio lobes of AGNs (Massa ro & Aiellol 
201 1), dark matt er decay (Olive & Silk 1985), and dark mat- 
ter annihilation (lAhn & Ko matsu 2005). It is not even clear 
that the observed background is real; it could be caused by de- 
tector backgrounds (e.g., caused by incident CRs or radioac- 
tive decays in the satellite). 

Not only the normalization but the spectral shape of the 
MeV background is mysterious. We plot energy flux per natu- 
ral log energy bin in Figure [Tjand similar figures; the area un- 
der a curve in an energy range in these figures is directly pro- 
portional to the power coming out within that energy range ^. 
It is clear that the MeV background represents more power 
than the GeV background. Furthermore, there are spectral 
breaks at ^ 3 MeV and ^ 20 MeV with unknown origin. 

We present here constraints and estimates of the star- 
forming galaxy contribution to the 7-ray background , with 
an emphasis on the MeV background. Soltan & Juchni ewic3 
( 1999) first suggested that starburst galaxies contribute to the 
MeV background through their IC emission, by scaling the 
IC emission from the observed synchrotron radio emission, 
assuming that starburst magnetic fields were the same as in 
die Milky Way (see also Soltan & Maixinkowski 2001). Most 
of the other estimates of the 7-ray background focus on pionic 
emission above 100 MeV; to the extent that leptonic emission 
is considered, the Milky Way spectrum is typically assumed. 
In the meantime, there have been a number of models of the 
7-ray emission from NGC 253, M82, and Arp 220 that take 
into a c count radio observations (iPaghone et al.lll996t iTorre: 
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MeV emission from star-forming galaxies and starbursts in 
particular are hampered by the very low sensitivity of current 
MeV instruments. Only weak constraints have been set on 
the MeV emiss ion from starburst galaxies with OSSE an d 
INTEGRAL ( B hattacharva"etan [19941; iDermer et al.|[T997[) . 
We model the MeV emission for starburst and normal galax- 
ies, and use the spectral shape to constrain the star-forming 
galaxy contribution to the MeV background. 

In Section[2] we review the MeV emission processes of star- 
forming galaxies. In Section [3] we present one-zone models 
of the Milky Way, M82, and NGC 253 that serve as templates 
for our 7-ray background calculation. In Section [4] we con- 
sider the effects of 77 cascades in the intergalactic medium 
on the shape of star-forming galaxy spectra. In Section [5] 
we present our results for the star-formation contribution to 
the MeV background, both the calculation of the MeV back- 
ground and constraints from the spectral shape. We assume 
Qm = 0.25, f^A = 0.75, and Ho = 70 km s"' Mpc"' . 



2. REVIEW OF MeV EMISSION PROCESSES 

Since 1//^ = E^dN/dE = EdN/d \nE, the area under the curve as plotted 
(X JuI^dlnE = J E^dN/dEdlnE = J EdN/dEdE. 
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Fig. 1 . — The observed unresolved extragalactic 7-ray background (blue) and the predicted 7-ray background from star-forming galaxies with our standard 
assumptions. Plotted data include the compilation in Gilh et al. (2007) (cross-hatching), SMM (shading), COMPTEL (X's), EGRET (triangles), and Fermi-hAT 
(squares). We show the total predicted background froin star-forming galaxies (black, thick solid), the total background from nonnal galaxies (black solid line, 
higher at MeV), and the total background from starbursts (black solid line, higher at 100 GeV). We also show various components of the 7-ray background, 
including the pionic emission (long-dashed) from starbursts and normal galaxies, positron annihilation radiation (long-dashed/short-dashed), nucleosynthetic 
7-ray lines and radiation from CCSNe (grey soUd), and cascade emission (from starbursts only; dotted). The separate contributions from IC and bremsstrahlung 
contiibution are not shown, but they make up most of the difference between the total and pionic emission. For comparison, the MeV background from type la 
supemovae from Horiuchi et al. ( 2011 ) is shown as grey shading. Our fiducial model for the Milky Way is the h = 2000 pc single zone model; for starbursts we 
use the M82 B = 300 /iG model; and we assume that 15% of the cosmic star-formation rate at all redshifts is in starbursts. However, the actual starburst fraction 
is highly uncertain and can affect the results at the order of magnitude level. 



At MeV energies, corresponding to temperatures > 10'" K, 
there are essentially no diffuse, thermal emission processes 
that can contribute to the 7-ray background. MeV emission 
from star-forming galaxies therefore either (1) involve non- 
thermal cosmic rays, or (2) directly comes from nuclear reac- 
tions, and sometimes both. We review the known MeV emis- 
sion processes that correlate with star-formation here. 

2.1. Leptonic CR Emission 

Most of the observed 7-ray emission from galaxies comes 
from CRs. The GeV-TeV emission is probably pionic, be- 
cause the injection rate of CR protons exceeds that of the pri- 
mary CR electrons by a factor ^ 50- 100 at GeV energies in 
the Milky Way. Even weak pionic losses as in the Milky Way 
are sufficien t to overwhelm the leptonic emission at these en- 
ergies (e.g.. IStrong et al.ll2010() . However, the kinematics of 
the pion production process leads to a rapid decline in the pi- 
onic luminosity below ~ 100 MeV (the "pion bump"), below 



which leptonic CR emission dominates. At a minimum, sec- 
ondary accompany pionic 7-rays; their leptonic emission 
can flatten out the pion bump ( Schlickeise jri982l) . 

Inverse Compton (IC) emission - Inverse Compton arises 
when a CR Compton scatters an ambient photon, boost- 
ing its energy in the observer frame. The IC emission is ex- 
pected to form a broad continuum, because the upscattered 
photon energy scales as the square of the electron energy (in 
the Thomson limit). It likely do minates the Milky W ay non- 
thermal continuum at 100 keV (iPorter et al.ll2008h . Below 
the Klein-Nishina cutoff, electrons of energy Eg will typically 
upscatter photons of energy e to AeE^ / {3meC^)^ . Ambient 
photons can come from the CMB, far-infrared (FIR) emis- 
sion from dust, or stellar UV/optical/infrared backgrounds; 
in the Milky Way, these all have similar energy densi- 
ties, while in starburst galaxies, the FIR emission domi- 
nates. Each of the radiation fields are roughly greybody fields 
with a temperature T and average photon energy ^ 2.7kT. 
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Therefore typical upscattered photon energies are w 
0.32(1 +z)£?o MeV, E^^ w 2A(Tpir/20 K)E^q MeV, and 
E°l' w 480 MeV(7;pt/4000 K)^^^ MeV for electrons of en- 
ergy lOEu) GeV upscattering the CMB, FIR, and optical ra- 
diation fields. IC from each of the components is expected to 
have a r = 7^/2+ 1/2 « 2.0-2.2 photon index, if the CR 
have a steady-state power law spectrum E~^ . 

The CMB energy density grows rapidly with redshift. IC 
upscattering off a strong high-z CMB can take energy away 
from other cooling processes, essentially quenching the radio 
luminosit y of a galaxy and transferring the power into MeV 
emission (ICaril li & Yun 1999; Murphy 2009). At high ener- 
gies, the CMB must overcome the Inverse Compton losses off 
the starlight and the synchrotron losses in the magnetic field. 
Thus, this effect is only important when Uqmb ^ Ub, or 

B<3.2nG(l+zf. (1) 
In Milky Way-like galaxies with B sa 6- 10 fiG and starlight 
energy density J/* « Ub, the CMB overwhelms the other 
losses before z 1. However, in starburst galaxies with 
B > 100 ^G, the CMB remains unimportant out to z > 5. At 
low electron energies, escape, bremsstrahlung, and ionization 
losses "buffer" the r adio spectrum against the C MB effect to 
even higher redshift (iLacki & Thompsonl2010ah . 

Finally, CR can escape into the halos of galaxies, where 
they cool by Inverse Compton emission. Indeed, Strong et al. 
(Ho 10) find that about - 30-60% of the power injected into 
CR electrons in the Milky Way is lost to Inverse Compton 
cooling, with greater IC luminosities for larger halos. The 
majority of this emission is below 100 MeV in most of their 
models, although the total < 100 MeV emission is still about 
a third or less of the total 7-ray emission dStrong et al.ll2010l) . 
Starburst galaxies are expected to be "electron calorimeters", 
with most o f the cooling in the starburst disk proper (e.g., 
|Volk'1989^, but radio haloes are observed in starbursts, in- 
dicating t hat some CR reach large heights off the plane 
dSeaquist & Odegardlll991l) . We do not consider these haloes 
further in this work. 

Bremsstrahlung - Bremsstrahlung radiation is emitted by 
CR when they are deflected by the electric field of a nu- 
cleus. A typic al bremsstrahlung photon has half the energy 
of the CR ( Schlickeiseill2002l) : thus the bremsstrahlung 
spectrum has the same shape as the CR spectrum, which 
is set by the dominant cooling or escape mechanism for CR 
e^. The bremsstrahlung loss time of an with energy 
Ecev GeV traversing pure hydrogen of number density n 
is fbiems ~ 40 Myr {n / cm~^)~\ but the ionization loss time 
is fion « 1380 Myr £'r,ev(w/cm"^)"'[ln£'r,ev+ 14.4]"' where 
Ec^y =E/GeV jStrong & Moskalenko|[l998h . Thus ioniza- 
tion losses therefore quickly softens both the CR and 
bremsstrahlung emission spectrum below ^ 400 MeV, in- 
dependent of density. Likewise, IC and synchrotron losses 
harden the CR spectrum at high energies, since they go 
as f (X £■"'. The result is a "bremsstrahlung bump" peak- 
ing at a few hundred MeV, where the energy-independent 
bremsstrahlung losses are most important. Thus, leptonic 
emission may itself be peaked near 100 MeV, even without 
hadronic emission. A low energy bremsstrahlung tail can be 
gene rated only if the injected elect ron spectrum is very steep 
(e.g.. lSacher & Schoenfeldeill983l) . 

Synchrotron and other emission from TeV-PeV e^ — Syn- 
chrotron emission from CR e^ deflected in a galaxy's mag- 
netic field is responsible for the low frequency radio emis- 
sion of galaxies, but it can also extend to X-rays and 7-rays. 



At MeV-emitting energies, there are few injected, but 
there are no competing cooling processes (IC is suppressed 
by Klein-Nishina effects). The typical energy of synchrotron 
photons is 1.2 MeV/ipgyBio, where Ep^y is the energy 
in PeV and Bio = 5/(10 /iG) is the magnetic field strength. 
In the Milky Way, the low magnetic field strength combined 
with the lack of IC emission observed at ^ PeV energies im- 
plies that there is little large scale hi gh-energy synchrotron 
emission (lAharonian & Atoyanll200dl) . The magnetic fields 
of starbursts are constrained to be much higher in starbursts 
like those in M82 and NGC 253 (-^ 100- 300 ^iG; e.g., 
[Thompson etani2006aL Ide Cea del Pozo et al. 2009b), so fliat 
synchrotron may be important at X-ray (iLacki & Thompson! 
2010b) and perhaps 7-ray energies. Synchrotron has a photon 
index of F > 2.0, steepening greatly at the tail of the CR 
spectrum. 

TeV CR e^ deflected in electrostatic waves will also 
emit MeV radiation through electrostatic bremsstrahlung, 
which o ccurs at higher frequencies for electrons of a given 
energy (ISchroder et al.l 120051) . The typical electrostatic 
bremsstrahlung photon energy for an electron of Lorentz fac- 
tor 7 is j^i^p compared to the synchrotron photon energy of 
-f^VB, where is the plasma frequency, i^b is the cyclotron 
frequency, and typically I'p ^ vb- However, electrostatic 
bremsstrahlung is only important if the energy density in elec- 
trostatic waves reaches the magnetic energy density; since this 
is not known, we ignore the process for the rest of the work. 

2.2. Positron Annihilation 

Nearly thermalized positrons can annihilate with ISM elec- 
trons into 511 keV photons, or they can form positronium 
atoms which can emit 3-photon annihilation continuum at en- 
ergies below 511 keV as well as 511 keV photons. A bright 
source of annihilation radiation is detected from the Galactic 
Center bulge, with a flux and spatial distribution not c onsis- 
tent with star formation (IWeidenspointner et al.l 1200^ ). but 
additional positron emission is detected from the Galactic disk 
('Weidenspo intner et al.ll20 08b'). The positron generation rate 
in the Galactic disk is consistent with a nucleosynthetic origin 
in th e fission of isotopes li ke ^^Al produced by massive stars 
(e.g.. iPrantzos et al]|201 lb . We consider this star-formation 
component of positrons, but not the (dominant in the Milky 
Way) bulge source of positrons in this paper 

In addition, CR positrons can annihilate while they are 
still relativistic (in-flight annihilation or lA), although most 
CR positrons survive to rest. Pionic positrons, produced 
alongside pionic 7-rays, a re injected at energies > 100 MeV 
(Beac om & Yuksell2006h . If we approximate the hadronic 
positron injection spectrum as QiE) = Q()5{E - Eo) where 
£■() ss 100 MeV, and since are cooled primarily by ioniza- 
tion (which has a cooling rate roughly independent of ener- 
gies) at £ < £0, the steady-state low-energy hadronic positron 
spectrum should have a nearly constant dN/dE. It can then be 
shown that the annihilation radiation fro m hadronic positr ons 
is hai-d, with dN-y/dE cx £'"'(ln(27)- 1) (lAharonianl2004 . so 
th at lyL^ ~ E^dN^/dE o c £'(ln(27)- 1) peaks near 100 MeV. 

IWatanabe et alJ (1 19991) calculated the 511 keV fine contri- 
bution to the 7-ray background, but did not include positron- 
ium continuum or inflight annihilation. 

2.3. Nuclear "f-Ray Lines 

Nucleosynthetic ISM j-ray lines - Some radioactive iso- 
topes synthesized in young, massive stars emit nuclear 7-rays 
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when they decay. The Al decay line at 1.809 MeV is the 
bright est 7-rav line above 1 MeV detected froni the G alactic 
plane ('Mahone v et al.|[T984t iDiehl et al.l [19951 l2006al) . The 
spatial distribution of the line is correlated with massive star- 
formation in the Galax y, includ ing individual star-form ing re- 
gions such as Cvgnus ([del Rio et al. 1996: Knodlse der et alJ 
I199I IDiehl et alJl2006air although older sources may con- 
tribute ^*A1 to the ISM. The only other detection of 7-ray 
lines above 1 MeV from the di ffuse Galaxy are the ^"Fe decay 
lines at 1.173 and 1.333 MeV dHarris et al.l l2005; Wan g etal] 
120071) . However, the 1.157 MeV '^'^Ti 7-ray line has been de- 
tected from the Cas A supernova remnant. The injection sites 
of Al and ^"Fe are thought to be supemo vae with some con- 
tribution from Wolf - Rayet star winds (e.g..|Knodlsede3ll999l : 
iPalacios et alJl2005l: lLimongi & Chieffi 2006). so their yields 
should scale with massive star-formation rate. However, some 
processes in old stellar populations, particularly n ovae, may 
also s ynthesize isotopes that emit 7-ray lines (e.g.. lJose et al.l 
[T997h . 

Gamma-ray lines from nuclear decay in the ISM are nar- 
row; Galactic observations of t he ^^Al line const rain line 
widths to be less than 3 keV jDiehl et al.ll2006bl) . but the 
line emission from a redshift distribution of galaxies will be 
smeared out into a continuum. 

The MeV background from these lines was also calculated 
in Watanabe et al. ( 1999). 

Core collapse supemovae - In addition to the long-lasting 
radioisotopes visible in the ISM, core collapse supemovae 
(CCSNe) themselves glow in 7-rays from short-lived iso- 
topes. The CCSNe are powered by the decay chain of 
^''Ni into 5^Co (half-life of 6 days), and then ^''Co into '^''Fe 
(half-life of 78 days). Each decay can produce 7-ray lines. 
However, CCSNe are highly Compton-thick, especially at 
early times, which downgrades much of the line emission 
into weaker continuum emission (e.g.. Ge hrelsetanfT987t 
lArnett & Fu 1989; Ruiz-Lapuente et al. 2001i). 

The 7-ray yield and MeV b a ckgro u nd from CCSNe was 
calculated by .W atanabe et al.l (Il999l). iRuiz-Lapuente et all 
(1200 Ih . andEwabuchi & Kumagail <200lh . 

Besides core collapse supemovae, there will be 7-ray 
line emission from type la supemovae, some of which are 
"prompt" after the formation of the progenitors. In this 
sense, some of the type la supernova emission can be as- 
sociated with star-formation. However, including this com- 
ponent would require a distribution of delay times of star- 
formation, which is beyond the scope of this paper In any 
case, the most recent predictions for this background show 
that it is only ^ 10 % of the observed MeV background 
dlwabuc hi & Kumag ail 120011; IStrigari et all 120051: lAhn et alj 
SosTHoriuchi & BeacomI MOm . We refer the reader to 
Horiuchi & Beacom (201(3) for discussion of the MeV back- 
ground from type la supemovae. 

CR de-excitation lines - CRs can also emit 7-ray lines 
through nuclear interactions with nuclei in th e ISM (e.g., 
iMeneguzzi & Reevesiri975l;lRamaty et alJI 19791) . In addition 
to narrow line components from the excitation of nuclei in 
the ISM, broader components are expected when CR nuclei 
themselves are excited. The strongest lines expected are the 
4.438 MeV '^C and 6.129 MeV "'O de-excitation lines (e.g., 
iRamatv et al.lll979l) . For these lines, the broad component 
should actually be more luminous than the narrow lines, be- 
cause CRs have much higher a bundances of C and O relative 
to H and He than the ISM (e.g.. lMever et alj|1998l) . 

CR-excited lines will be much fainter than the pionic 



emission unless there are many low energy CR nuclei, 
with only 0.01-0.1 line ^ MeV photon but ^ 1 pionic 
phot on per ~ 100 M eV/» nucleus if losses are unimpor- 
tant dParizot & Leho ucq 2002). A low energy component 
of CR nuclei, peaking at only a fev v MeV, is sugg ested 
by some CR ionization studies (e.g., Ilndriolo et al.l 12009,) . 
but CR- excited nuclear lines have no t been detected in the 
Galaxy dTeegarden & Watanabell2006l) . A claimed detection 
of the 4.438 M eV and 6.129 MeV l ines in the Orion star- 
forming region ([Bloemen et al.lll994l) later proved spurious 
dBloemenetal.ll 19991) . 

2.4. Emission from Discrete Sources 

Star-forming galaxies are unresolved spatially in 7-rays, ex- 
cept for the Milky Way and the Magellanic Clouds. How- 
ever, in addition to the diffuse 7-ray emission, compact 7-ray 
sources as observed in the Milky Way may contribute to the 
total galactic 7-ray luminosity. Many of these sources - in- 
cluding supernova remnants, pulsar wind nebulae, and molec- 
ular clouds located near CR accelerators - are directly related 
to star-formation, and could increase the contribution of star- 
forming galaxies to the 7-ray background. A study by IN- 
TEGRAL indicates that at > 1 MeV, most of the Galactic 
emission is an unresolved diffuse compo nent, unlike at lowe r 
energies where point sources dominate (Bouchet et al. 2008). 
However, given how poorly the MeV regime is understood 
even in the Milky Way, with even the advanced GALPROP 
models failing to account f or all emission in the ^ 30 MeV 
range (Strong et al. 2004bl iPorter et al.]|2008b . it is possible 
that additional MeV emission will come from unresolved dis- 
crete sources (c.f., Strong_et al. 2000b). We do not consider 
discrete sources further here, except to note them as a great 
uncertainty. 

2.5. Brief Summary 

Over much of the range a few hundred keV to a few hun- 
dred MeV, IC is expected to be the dominant emission pro- 
cess. This is because it forms a broad continuum; further- 
more, both radiation fields and of the relevant energies 
(GeV-TeV) are known to be present in star-forming galaxies. 
Nonthermal bremsstrahlung may be important at ^ 100 MeV, 
but it appears as a bump in the spectmm that subsides at lower 
energies. Line emission from CCSNe, positron annihilation, 
and ^^Al may be important from a few hundred keV to an 
MeV. The tail of the synchrotron spectmm also may be im- 
portant at a few hundred keV, though its strength is highly 
uncertain. 

3. MODELLING OF 7-RAY SPECTRA 

To estimate the MeV - GeV background contribution from 
normal and starburst galaxies, we need to understand the 7- 
ray SED of star-forming galaxies - both the spectral shape 
and the normalization per unit star-formation. Since much of 
this emission is from CRs, this requires an understanding of 
CR populations in both normal and starburst galaxies. The CR 
population is governed by the diffusion-loss equation, which 
includes terms for CR injection, spatial propagation, cooling, 
and destmction (see Appendix). However, for the purposes of 
understanding other galaxies, we can reduce the complicated 
diffusion-loss equation into a one-zone leaky-box equation, 
which treats the galaxy (or its starburst core) as one region 
with an homogeneous environment and CR population. 

There are very detailed models of the nonthermal broad- 
band spectmm of the Milky Way, particularly with the 



6 



LACKI, HORIUCHI, & BEACOM 



GALPROP code (e.g.. IStrong et alj|2010l) . and several one- 
zone and 3D models of the st arburst galaxies M82 and NGC 
253 dPaglione et alj Il996t Domingo-Santamaria & Torres' 
20051: Persic e t al.l l2008t Fd e Cea del Pozo et al. 2009a; 
Rephaeli et al. 200^. iLacki et al.l 1201 If) . However, these 
models use the present CMB; at high z, the strong CMB can 
increase Inverse Compton losses and change the broadband 
spectrum. In addition, these models generally do not include 
positron annihilation or nuclear lines. We therefore run our 
own models, using the Milky Way as a spectral template for 
normal galaxies, and M82 or the NGC 253 starburst as a 
spectral template for starbursts. 

3.1. Procedure for One-Zone Models 

We solve the steady-state, one zone CR propagation equa- 
ti on (l eaky box model) using a Green's function described 
in [Torres (2004). which we solve with a code described in 
ILacki et al.] jloWtt . We include primary nuclei and e*, sec- 
ondary (knock-on and pionic) e^, and higher order pair 
from 77 annihilation. Our procedure is largely the same as in 
iLacki & ThompsoU (l2010 b). although with several new fea- 
tures: we include positron annihilation, nuclei heavier than 
hydrogen (both in the primary CRs and the ISM), and nuclear 
de-excitation lines (as described in the Appendix). 

We model star-forming galaxies as disks with radius R and 
midplane-to-edge scale heights h. The CR injection power is 
proportional to the supernova rate: a fraction 77 ^ 0.1 of the 
10'''£'5i ergs in supernova kinetic energy goes into primary 
protons, while a fraction ^ ^ 0.01 goes into primary electrons. 
For the starbursts, we relate the supernova rate directly to the 
TIR (total infrared: 8- 1000 jim) luminosity, which should 
scale almost directly with star-formation rate: 



r^S= 0.036 yr-'Vi7 



^TIR 



10 



10.5 



(2) 



where tAiv = 1 is an initial mass function (IMF) de- 
pendent constant ([Thompson. Ouataert. & WaxmanI 120071: 
iLacki (fe ThompsoU 1201 Obi) ! Since massive stars domi- 
nate both the bolometric luminosity and the supernova 
rate, the supernova rate inferred from the luminosity 
should vary by only a few percent for different IMFs 
([Thompson. Ouataert. & Waxman|[2007h . For the Milky Way, 
we relate the supernova rate to an assumed star-formation 
rate: 

SFR 



- SN 



: 0.0084 yr-'esaiA^-iv 



M. 



yr 



(3) 



where esaiA = e/(4.9 X 10 "*) is a conversion between bolomet- 
ric lumin osity and star-formation rate corrected to a "Salp eter 
A" IMF ( Kennicutll [199^ : [Baldry & G lazebrook 20 031), as 
used for the cosmic star-formation rate in Hopkins & Beacorrj 
(^006).^ From the supernova rate, we get the volumetric 
power injection for CR protons: 



ecR,p = 10-' 'ergs EsiriT^nl{2nR^h) 
and primary e^: 

ecR,e = lO^'ergs Esi£V^^/{2T:R^h) 



(4) 



(5) 



' IHopkins & BeaconJ f200^ find that for a given bolometric luminosity, 
a "Salpeter A" IMF has only 77% of th e star-formation rate of a standard 
Salpeter IMF from 0. 1 - 100 Mq used in lKennicutll fT993) . 



Primary CRs are injected with a dQ/dq = Cq~'' power 
law spectrum, where q is momentum: this is the test- 
particle approximation, although nonlinear physics in ac- 
celeration may alter the inject i on spectrum at low en- 
ergies (e.g., Berezhko & ElHson 1999; Elli son et al.1 [20001: 
[Malkov & O'C Drury 2001: Blasi et al...2005i) . The normal- 
ization for protons and e^ is set by: 



ecR = C / q ''(\/ q^c^ + irP-c'^-mc^)dq, 



(6) 



where m is the rest mass of the accelerated particle, ^^in and 
q'max are the momenta corresponding to the minimum and 
maximum kinetic energy particles are injected with, and ecR 
is the volumetric power injected in CRs of a given species. 
We use a minimum injection kinetic energy of 1 MeV; for 
the protons and nuclei, we use a maximum injection Lorentz 
factor of lO*" (K^^^ « 1 PeV), while for e±, the maximum m- 
jection Lorentz factor is either 2 x 10^ (^^nax ~ 1 TeV) for the 
Milky Way, or a free parameter for the starbursts. With these 
injection spectra, the ratio of dQ/dq for primary protons and 
electrons at energies > nipC^ is approximately: 



Krr, 



p-1 



(7) 



which i s known to be ~ 5 0-100 in the Milky Way (e.g., 
Ginzb urg & Ptusk in 1976j [Warren et all [20051: Istrong et al.l 
i2010i) . The normalization for the nuclei are scaled using the 
Milky Way CR abundance ratios (see the Appendix). 

We run models for hadronic-origin CRs (nuclei, pionic sec- 
ondaries, and the radiation and 77 pair they generate) 
and leptonic-origin CRs (primary CR and the radiation 
and 77 pair they generate) separately, and scale them by 
77 and ^ respectively. These components are then added to- 
gether Diffusive losses are applied to CRs of all species 
in both the Milky Way and starbursts. Advective losses are 
also applied in starbursts. CR nuclei experience catastrophic 
inelastic (including pionic) losses and continuous ionization 
cooling. CR electrons experience ionization, bremsstrahlung, 
synchrotron, and IC cooling; CR positrons experience these 
losses as well as catastrophic annihilation losses. We include 
secondary from knock-off processes and pion production 
in nucleus-nucleus collisions, as well as tertiary e^ from 77 
pair production. Secondary nuclei from collisions, however, 
are not included. We calculate the emission from pionic 
collisions (both 7-rays and neutrinos), synchrotron, positron 
annihilation, several of the strongest nuclear lines, IC, and 
bremsstrahlung. More technical details are given in the Ap- 
pendix. 

The propagated spectra depend on a number of param- 
eters, many of which are poorly constrained in starburst 
galaxies. The GeV-TeV 7-ray emission, which is pionic in 
most of our models, constrains the proton acceleration ef- 
ficiency r]. The radio data constrains a combination of the 
electron acceleration efficiency ^ and B, but there is a de- 
genera cy between these paramete rs (see t he similar discus- 
sion in | Lacki & T hompson'2010b'; also see 'Persic et al. 20081 
[de Cea del Pozo et al. 2009 a. and RephaeU et al. 2009) . With 
high B, few e^ are needed to produce the observed GHz syn- 
chrotron radio emission, and the population in starbursts 
is dominated by hadronic secondaries. Models with Milky 
Way values for the proton/electron injection ratio have high B. 
Models with low B have large ^ and therefore large amounts 
of leptonic 7-ray emission in the MeV bands. In addition, the 
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the CR spectrum. Finally, we add a component of free-free 
emission with spectrum Si^ = 5ff(l GHz)(j//l GHz)"" ', which 
flattens out the radio spectrum. 

To keep the number of spectral templates for the starbursts 
to a manageable level, we group models by B. We consider 
consider magnetic field strengths of 50, 100, 150, 200, 250, 
300, and 400 for M82 and NGC 253's starburst. For each 
B, we minimize for the sum of radio and 7-ray data for the 

par ameters 77, Vwind, h 7S™, %(! GHz), and fdiff(3 GeV) 
(c.f. lCrocker et al.l201 Ibl) to select a best model for that B, for 
both primary electrons and an equal mix of primary electrons 
and positrons (listed in Table [TJ. In the Milky Way, most of 
the parameters are constrained by direct observations and we 
do not do fitting. 

To the total 7-ray spectra, we then add in the nucle- 
osynthetic 7-ray lines. The 1.809 MeV line from -^Al is 
scaled from its Milky Way luminosity using the supernova 
rate: Li.go9 = 1.7 X 10"ergs s-i(rsN/rsN,Mw) (Di ehl et alJ 
I2006ah . In the Milky Way, the 1.173 and 1.333 MeV 
^''Fe lines are observed to have 0. 11 photon each fo r ev- 
ery photon of the 1.809 ^^Al fines dHarris et al.ll2005b . We 
therefore set Li 173 = 0.11(1. 173/1.809)Li §09 and Li 333 = 
0.11(1.333/1.809)Li.809- Finally, we scale the 1.157 MeV 
^"^Ti 7-ray line directly with the supernova rate: L 1.1 57 = 
1.6 X 10^**(M7-,_44/10-'' Mq) ergs s-i(rsN/yr"'), where M7-,-44 
is the yield of ^''Ti per supernova. The '*^T i yield is controver- 
sial, bec ause Cas A is br ight in the line (llvudin et al.lll994t 
see also IVink et al.ll200TI for detections of associated X-ray 
lines) and S N 1987A's emission indicates large yields (e.g.. 



wind speed Vwind, gas scale height h, primary electron cutoff Lqmi = 5.9 x 10 ergs s (FsN/rsN.Mw) and L1.238 = 6.7 x 
7m™, and diffusive escape time fdiff(3 GeV) can all influence lO-'^ergs s"'(rsN/rsN,Mw) (iRuiz-Lapuente et alJ 1200 Ih . In 

addition, we include the continuum emission using the spec- 
trum calculated in Ruiz-Lapuente et al.. (.2001,) . scaled to the 
0.07 M0 of ""^Ni per supernova. 

3.2. Milky Way 

Introduction - The Milky Way (MW) is the only non-active 
galaxy observed in MeV 7-rays, and only one of a handful of 
normal (non-starburst, non-active) galaxies observed in GeV 
7-rays (with the LMC, SMC, and M31). At low energies, 
below ~ 50 keV, most of the Galactic X-ray ridge emission 
is resolved into stellar s ources, such as cataclysmic variables 
jRevnivtsev et al.|[2009h . The energy region up to 511 keV 
includes a spectral spike from positron annihilation radiation. 
This annihilation radiation is concentrated towards the cen- 
ter of the Galaxy, which is not expec ted if it is correlate d 
with star-formation (see the review by iPrantzos et al.ll201 lb . 
In the energy range from 100 keV to 100 MeV, there ap- 
pears to be a power law continuum, with F w 1.5, thus in- 
creasing in luminosity to higher energies. This power law 
continuum has been interpreted in the pa st as bremsstrahlung 
(ISchlickeiseiill982t ISacher & Schoenfeld er 1983) , though it 
is now believed to be Inve rse Compton emission ( Porter et al.l 
l2008tlBouchet et al]|201 lb . However, there also appears to be 
an excess of 7-rays above the power law detected by COMP- 
TEL in the 7^ 5 - 20 MeV range towa rds the inner Galaxy 
(ISti-ong et al.ll2004bl: iPorter et al.ll2008l) . The Galactic 7-ray 
luminosity peaks around a GeV, a pionic feature e xpected 
from t heoretical exp ectati ons of CR propa gation (e.g. JSteckerl 
fT970t lHunter"etaLlfT997t IStrong et al.l l2004b). Finally, the 
emission drops off rapidly as E~^-^ to E~^-^, continuing to 
TeV energies in at least some regions of the Galactic Plane 
(iProdanovic et al.ll2007HAbdo et al.ll2008h . 

In modeling the Milky Way, we use a radius of 10 kpc. We 
consider scale heights (midplane to edge) of 1 kpc and 2 kpc. 
The gas surface density is measured to be ^ 0.0025 g cm"^ 
in the Solar neighborhood and through much of the Galactic 
disk (.Boulares & Coxlll990l: lYin et al.ll2009l) . For the radia- 
tion field, we average the maxirn um metallicity gradient radi- 
ation fields of lPorter et all (|2008|) for 7? = 4 kpc and = 8 kpc, 
reflecting the fact that much of the star-form ation (and there - 
fore CR generation) occurs at ss 5 - 7 kpc (lYin et al.ll2009h . 
In addition, we use a magnetic field strength of 6 /iG (e.g., 
IStrong etaTI l2000at iBeckl l20l"Tl) . a diffusive escape time of 
friiff(£ ) = 30 Myr (E/3 GeV)-!/^ (IConnellll 199 8: Webbe r et all 
lIOOl . and a maximum primary electron Lorentz factor of 
7m™ = 2 X 10^ as inferred from the apparent cutoff in the pri- 
mary electron spectrum at TeV energies (e.g.. lAharonian et al.l 
2008, 2009). 

We assume that the total Galactic star-formation rate is 
2 Mq yr~', giving us a supernova rate of 1.6 century"' for 
the Sal A IMF. Estimates of the star-formation rate and su- 
pernova rate of the Galaxy span nearly an order of magni- 
tude, though most of the larger ones are from older works 
(see Table 1 of iDiehl et alj l2006ah . Recent estimates 
of the star-formation rate from free-free emission or di- 
rect counting of young stell ar objects find star-form ation 
rates of 0.9 - 2.2 Mp^ yr"' dM urrav & Rah maiil 120101) and 
0.7-1.5 Mq yr"' (iRobitaille & Whitney 20 Irrespectively; 
therefore, our star-forma tion r ate estimate may be some- 
what high. Chomiuk & PovichI (|201 1.) compile several star- 
formation rate estimators for the Galaxy and conclude that it 



Chugai et al. 1997; Fransson & Kozma 2002; Jerk strand et al.l 
20111). but~the Galaxy as a whole is faint in 1.157 MeV 
(iThe etal.ll2006h . We assume 10 ^ Mq per supernova, but 
estimated values range from 10"^ -2 x 10"^ M©. 

In addition, we run a one-zone model of nucleosynthetic 
e"^ from radioactive isotopes injected as a delta function at 
1 MeV, scale d to the positron gener ation rate in the Galac- 
tic disk from IWeidenspointner et all (|2008b): (dQ/dE)s^ = 
8 X 10'*2s"'(rsN/rsN,Mw)5(^^-l MeY)/ (InR^h). The 7-ray 
spectrum is then the sum of the emission from hadronic-origin 
CRs, the leptonic-origin CRs, the nucleosynthetic e"*", and the 
nucleosynthetic 7-ray lines. 

When we have a best-fit model for each B and primary 
composition, we then use those same parameters run models 
at redshifts from 0.0 to 4.9 at intervals Az = 0.1. The only 
alteration with redshift is that the radiation field is altered so 
that the CMB is increased to its value at each redshift. No 
other evolution of the spectrum shape of each class of galaxies 
(due to, e.g., increased gas fraction at high z) is considered. 
Instead, galaxy evolution is assumed only to convert galaxies 
from being starbursts at high z to normal at low z (see the 
extensive discussion in section 15.2b . besides the increase in 
cosmic star-formation rate of both kinds of galaxies at high z- 

The 7-rays from CCSNe are not added to the mod- 
els of the individual galaxies (Milky Way, M82, and 
NGC 253's starburst), since supemovae are rare and tran- 
sient events. However, we do include them in the 7-ray 
backgrounds. We assume in our calculations tha t each 
CCSN produces 0.07 M^ of ^'^Ni (e.g., Catcht^o le et all 



1988t lArnett&Fulfl989l: iBouchet & Danziger»1993i: iHamuyl 
2003h . We include the strongest 7-ray lines from ^^Co at 



0.847 MeV and 1.238 MeV, with time-averaged luminosities 
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Fig. 2. — The predicted 7-ray luminosity spectrum of the Milky Way (solid), compared to extrapolated observations of the inner Galaxy with SPI (pentagons), 
COMPTEL (squares), and EGRET (triangles), as well as GalProp models (grey shading; Strong et al. 2010). The different components to the emission are 
positron annihilation (long-dashed/shoit-dashed), nucleosynthetic 7-ray lines (grey solid), CR nuclear de-excitation lines (grey solid), bremsstrahlung (short- 
dashed), IC (dotted), and pionic emission (long-dashed). The excess of GeV emission is an instrumental artifact of EGRET. Note that the predicted positron 
annihilation rate only includes the contribution from star-formation, but the observed positron emission includes the Galactic Center source. Emission from 
CCSNe is not included. Also note that the observations, predicted line spectra, and total predicted spectra have different energy binnings; the binning for the total 
spectrum is C^EjE = 3%. 



is ~ 1 .9 M0 yr"' for a Kroupa IMF. For a given star-formation 
rate, and for stellar masses greater than 1 where most of 
the starlight used to infer the star-formation rate history of the 
Universe comes from, the IMFs differ by less than 30%. At 
10 Mq, the approximate mass of most supernova progenitors, 
the IMFs differ by only 13%. The result of 2 M© yr"' should 
thus be applicable for the SalA IMF that we use, to within a 
factor of < 2. On the other hand, we use the star-formation 
rate solely to calculate t he supernova rate, w hich is about a 
factor ^ 1 .2 lower than in lDiehl et al.l(l2006ab (estimated from 
the ^''Al 7-ray line). 

We set the 1 MeV positron injection rate from radioactive 
elements for the entire Galax y to be 8.1 x lO'^^s"^ as deter- 
mined for the Galactic disk byl Weidenspointner et al.l (l2008ah 
with their "bulge -1- thick disk" model. This does not include 
the positron rate from the bulge, which is substantially higher 
(~ 2 X lO^^s"'). In using this value for the nucleosynthetic 
positrons, we assume that CR positrons contribute a minority 
of this r ate, which is verifie d by our modeling and GALPROP 
models jPorter et al.ll2008l) . 

Comparison with observed ^-rays - We scale our predicted 
7-ray spectra by comparing our models' 100 MeV - 100 GeV 
pionic 7-ray lumin osities to the more detailed predictions in 
IStrong et all dMTol) . We find using standard acceleration ef- 
ficiencies (77 = 0.1) and our fiducial supernova rate that our 
models' 7-ray luminosity is too small. We find we must scale 
up r] (and ^ using 5 = 15) the h = 1000 pc 7-ray luminosity by 
a factor 1.7 and the h = 2000 pc luminosity by a factor 3.3. 
There could be several reasons for this discrepancy: our su- 
pernova rate could be too small by a factor ^2-3 (but this 
would contradict current observations), the CR acceleration 
efficiency could be greater by a factor ^ 2 - 3, or the strength 
of pionic losses to escape losses for CR nuclei is ~ 2 - 3 times 



stronger than we calculate (meaning either slower escape or 
higher gas densities). 

To check these reasons for the low 7-ray luminosity, we ex- 
amine the density and CR luminosity of our models. The aver- 
age density that CRs experience in our models is 0.24 cm~^ in 
the h = 1000 pc model and 0.12 cm~^ in the h = 2000 pc model; 
CRs observed at Earth are inferre d to traverse mean densities 
of 0.2 5 cm-3 in the Milky Way (IConnel]|ll99§: ISchlickeised 
I2OO2I) . The h = 2000 pc model therefore is too 7-ray dim 
partly because it is too low density; this is not the case for 
the h = 1000 pc ca se. The CR l umino sity in our models is 
5.1 X 10"*° ergs s"': IStrong etaP (l2010l) find CR luminosities 
of 7 X 10 ergs s"' for their h = 2 kpc models. Since we 
are comparing to their 7-ray luminosity, our CR luminosity 
is ^ 40% too low; after accounting for this, the 7-ray lumi- 
nosity of the h = 10 00 pc model is within 25% of that of the 
IStrong et al.l (1201 Ol) model. We conclude that the need for the 
scale factor arises because the CR luminosity is too low in our 
models, and also because the mean density is too low in the 
h = 2000 pc model. 

The rescaled 7-ray spectra are plotted in Figure |2] Pionic 
emission dominates the 7-ray emission at most energies above 
100 MeV, leading to the observed pionic bump. At MeV ener- 
gies, the relative strengths of the leptonic emission processes 
depend on the assumed scale height. In the smaller scale 
height (1000 pc; left panel in Figure|2]l models, CR expe- 
rience higher gas densities than the large scale height models 
(2000 pc; right panel in Figure |2]i, but the radiation field is 
the same. Therefore, bremsstrahlung receives more power in 
the h = 1000 pc model than the h = 2000 pc model. We find 
that in the h = 1000 pc models, bremsstrahlung is the most 
important emission from ^ 10-100 MeV and IC is the most 
important continuum emission at energies below that. In the 
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h = 2000 pc models, IC is the dominant continuum emission 
process at energies below 100 MeV, although bremsstrahlung 
nearly equals it at 100 MeV. This is in accordance with ex- 
pectations that the truly diffuse Galactic nonther mal contin- 
uum in the 100 ke V - MeV range is IC emission ( Porter et al.l 
I2OO8: Bouchet et alJl20lTl) . The /z = 1000 pc model spectrum 
is essentially in agreement with GALPROP luminosity spec- 
tra within uncertainties, although the h = 2000 pc model h as 
much more leptonic emission (Figure l2l lsFong et al.ll20Toh . 

In Figure |2] we also plot observations of the inner Milky 
Way 7-ray s pectrum. The shown data are scaled from INTE- 
GRAL SPI feouchet et alJ |200i) COMPTEL and EGRET 
(iPorter et al.ll200 8') data plotted inl fUrler et al.l (l2010l) . These 
are translated into Galactic luminosities by assuming that CR 
7-ray emission is distributed spatially the same as ^^Al nu- 
clear line emission (which is associated with massive stars), 
and using the ratio of the ^^Al line flux observed from the 
same sky region a nd the Galactic luminosity reported by 
iDiehl et aU (l2006ah . Comparing the spectral shape of the 
models with the observed spectrum, we find that the li = 
1000 pc model is the better match below 1 MeV (see Fig- 
ure |2]l. The IC continuum in the h = 2000 pc model is up 
to a factor ~ 3 times too high at these energies. Conversely, 
the h = 2000 pc model is the better match above 1 MeV, with 
the h = 1000 pc model underproducing flux in the few MeV 
rang e by up to a facto r ^ 3. Note the GALPROP models 
from lStrong et al.ll2010l also fail to reproduce the MeV emis- 
sion, even considering the uncertainties in halo size. The 
plotted data above the GeV are from EGRET, and the appar- 
ent "GeV e xcess" is now known to be an instrumental effect 
dStecker et al. 2008; Abdo et al. 2009). 

We use the h = 2000 pc model as our fiducial model of the 
Milky Way. 

Results for MeV j-rays - Most of the star-formation asso- 
ciated positron annihilation rate comes from nucleosynthetic 
positrons. Of the 8.1 x 10^^ s"' nucleosynthetic positrons in- 
jected in the Galactic disk at 1 MeV, we find that 6.8 x lO"*- s"' 
(83%) survive to low energies. In addition, we find that 
7.6-7.8 X 10^' s"' positrons are injected through pion pro- 
duction, about 40-80% of the rate calculated by GALPROP 
models (Porter et al. 2008). Of these, 5.5 x 10"*' s"' (72%) an- 
nihilate near rest in the h = 1000 pc model and 4.8 x 10^' s"' 
(61%) annihilate near rest in the h = 2000 pc model. The re- 
maining positrons either escape or annihilate in flight. Our 
models indicate that ^ 7% of the disk positrons are hadronic 
in origin. The luminosity in positronium continuum and the 
511 keV line is ^ 6 x lO-'^ ergs s"', whereas the i^L^ luminos- 
ity of the nonthermal Galaxy at that energy is 4 x lO^*^ ergs s~' 
in the h = 1000 pc model and 7 x 10 ergs s"' in the h = 
2000 pc. Therefore positron annihilation near rest associated 
with star-formation is only ^ 10% of the luminosity of the 
Galaxy near 511 ke V. 

The total positron in-flight annihilation luminosity is one 
half of the luminosity of positrons annihilating at rest, ^ 3 - 
4 X 10"'^ ergs s~'. While by number, most positrons survive 
to low energy, pionic positrons have much more energy when 
they are injected than when they are near rest. Therefore, most 
of the power from pionic positron annihilation comes out near 
^ 65 MeV. On the other hand, the other nonthermal Galactic 
emission is also greater at these energies, with i^L^ w 1 .2 - 
1.8 X 10^^ ergs s"' so that in-flight annihilation is just a few 
percent of the Galactic luminosity at these energies. 

The nucleosynthetic 7-ray lines are a significant but minor- 
ity contributor to the MeV 7-ray luminosity of the Galaxy. 




Fig. 3. — Rest-frame evolution of the total (black) and IC (grey) 7-ray 
spectrum of the Milky Way (with its cun'ent physical conditions) placed at 
different redshifts (0: solid; 1: long-dashed; 2: short-dashed; 4: dotted). At 
z = 0, the emission falls rapidly as decreases. But at higher redshift, the 
stronger CMB enhances Inverse Compton emission in the MeV band until it 
is at least ~ 1/2 of the GeV luminosity, paitly flattening the pion bump. To 
show the differences in spectral shape, we assume the Milky Way's SFR is 
the same at all z, although it was likely bigger in the past. 

The strongest line, 1.809 MeV emission from ^*A1, has a lu- 
minosity of 1.7 X 10^^ ergs s"'. The other nucleosynthetic 
fines considered (''"Fe 1 . 173 and 1.333 MeV, '^'^Ti 1 . 157 MeV) 
have luminosities of 1.5-2.5 x 10^^ ergs s~', for a total nu- 
cleosynthetic line luminosity of 2.6 x 10^^ ergs s"'. The non- 
thermal luminosity at 1.5 MeV is i/L^ 



4 X 10 ergs s"' for 

h = 1 kpc and 8 x 10^*' ergs s"' for h = 2 kpc. Nucleosyn- 
thetic lines therefore make up about ^ 30% of the luminosity 
at energies around 1 MeV. 

More exotic processes, such as synchrotron and nuclear de- 
excitation lines, are insignificant in our models of the Milky 
Way. 

High z - In Figure [3] we see the effects of the increasing 
CMB strength at high z: the Inverse Compton component be- 
comes much stronger This enhances the MeV emission con- 
siderably in Milky Way-like galaxies: in fact, at z = 2, the 1 
MeV I'Liy luminosity is 76% of the 1 GeV i^L^ luminosity. 
Thus the enhancement of IC at high z must be accounted for 
in calculating the MeV background from MW like galaxies. 

An interesting effect that occurs in Milky Way-like galaxies 
at high z is that the fraction of hadronic positrons surviving to 
rest increases from 60-70% to ^ 87%, enhancing the hadronic 
positron annihilation rate by ^ 20%. This is because the in- 
creased Inverse Compton losses off the CMB grow stronger 
than escape, so that CR positrons are trapped more effectively 
in the host galaxy. However, given that hadronic positrons 
form only a small fraction of the star-formation associated 
positrons, the total Galactic disk positron annihilation rate is 
only enhanced by 2%. 

3.3. M82 

Introd uction - M82 is a n earby starburst galaxy {D « 
3.6 Mpc: lFreedman et al.ll 19941) . the brightest in flie GeV and 



10 



LACKI, HORIUCHI, & BEACOM 



TeV 7-ray sky (lAbdo et all 1201 Oft lAcciari et all l2009h. and 
one o f the brightest galaxies observed in the IR dSanders et alJ 
l2003h . The starburst lies in the central regions of the 
galaxy (R < 250 pc), where a large amount of gas, infrared 
emission, and radio emission is observed ( Goetz et al.lll990l: 
IWilUams & Boweiil2010) . About 2 x 10** M© of gas is ob- 
served in the c entral regions, giv ing a gas surface density 
of 0.17 g cm-2 dWeiB et al.ll2"00lb . M82 is also observed to 
host a starburst wind, and we consider wind speeds from to 
lOOOkms-i. 

Usi ng the observed TI R luminosity of 5.9 x 10'° L© 
from (ISanders et al.l l2003h . we take the supernova rate as 
0.067 yr"', compared to typically quoted values in the liter- 
ature of ~ .1 yr~' with facto r ^2-3 uncertainty (see the 
summary in Lacki e t al.ll20l"ll) . This corresponds to a star- 
formation rate of 8.0 M© yr"'. 

Essentially nothing is known observationally about M82's 
MeV emission. It has been detecte d in 10-50 keV X-ray 
emission by Suzaku and Swift-BAT (ICusumano et alJl2010l) . 
Very weak upper limits on the continuum MeV emission ex- 
ist from non-detections by OSSE (Bhattacharya et al. 1994). 
Searches for nuclear lines from M82 and NGC 253 with 
OSSE have also found nothing yet (iBhattacharva et al.lll994l) . 
In the two year Ferm/-LAT catalog, M82 is detected with 
> 5(7 significance at energies down to 100 MeV, although the 
ferm/-LAT data does not yet go down below that. There are, 
however, several theoretical models for M82's multiwave- 
length nont hermal e mission, including MeV e nergies (e.g., 
iPersic et al. 2008; d e Cea del Pozo et al.ll2009ah . 

Constraints from radio and GeV-TeV ^-rays - We use 
the int erferometric radio data compiled in William s & Bower! 
(1201 Ol). combined wi th the original 7 -ray data from 
lAbdo et al.1 (l2010fl) and lAcciari et all ( 120091) . The combined 
radio (Fig. |4]l and 7-ray (Fig. |5]l data constrain the allowed 
parameter space in M82 significantly, although there are still 
large degeneracies. We also plot X-ray data on Figure |5] 
although we do not include it in the fitting. This emission 
largely comes from discrete X-ray binaries, and even the dif- 
fuse emission likely comes from hot t hermal gas which would 
not co ntribute at MeV energies (e.g., [Strickland & HeckmanI 
l2007h . 

The best-fitting Z? = 50 //G model (which has h= 100 pc) 
is a relatively poor fit to the data compared to models with 
higher B (the total for B = 50 fiG models is > 120, whereas 
> 70 for models with higher B). This is because the ob- 
served radio emission combined with the assumption of a 
small magnetic field imply a large number of CR e^, in turn 
implying vast amounts of leptonic 7-ray emission which ex- 
ceed the observed amounts. Although there may be other pa- 
rameters affecting the emission, the likely explanation is that 
B is greater than 50 ^iG in M82. 

The best-fit models for other B (100-400 ^lG) all fit M82 
with comparable accuracy (x^ = 70-73) to each other The 
best-fit models all have relatively small scale heights of 30 
pc, or in the case of B = 400 ^G, 10 pc (see Table [T]i. This 
is because of the shallowness of the radio spectrum com- 
bined with the spectral curvature, which is inconsistent with 
large amounts of free-free emission ( Williams & Bower 2010.: 
Fig. IDl. Smaller h for a given Y,g implies higher gas densities 
and therefore stronger bremsstrahlung and ionization losses; 
the smaller h also means the wind-crossi ng times are shorter 
for M82 (compare the expectations in iLacki & Thompson! 
llOlOb,) . In all of the B > 50 ^lG models, %(1 GHz) = 
0.34-0.45 Jy for a thermal fraction of ^ 0.04 (compared to 



100 



M82 

B = 300 /j,G 




100 



1/ [GHz] 



Fig. 4. — Predicted radio spectrum of M82 in the best-fit B = 300 fj-G 
model witli no primary positrons. Plotted lines are total radio emission 
(solid), total synchrotron emission (long-dashed), synchrotron from secon- 
daries (short-dashed), synchrotron from primaries (dotted), and free-free 
emission (dash-dotted). Plotted data are the interferometric data compiled 
in Williams & Bower (2010) (open circles), and data from Klein et al. 1 198^ 
(X's). The bottom panel shows the residuals to the radio data. The discrep- 
ancy at low frequencies (data at < 1 GHz is not included in our fitting) is 
probably due to free-free absorption. At high frequencies, thermal free-free 
emission flattens the spectrum. 



- 0.06 in !Wimams & Bowei!l20Tol) . 

The relatively poor 7-ray data do not allow us to discrim- 
inate the relative mixture of leptonic and hadronic emission 
at GeV energies. The GeV emission in B = 100- 150 /iG 
models is actually dominated by leptonic emission, both 
bremsstrahlung and Inverse Compton. These models have 
very high ^ (0.03-0.05) and small r] (0.009 - 0.03), resulting 
in high (5 w 1 . However, in order for the flat Inverse Compton 
continuum to not overproduce the TeV emission, the primary 

electron spectrum must cut off at 7m"x = 10'' to 10^. The TeV 
emission is pionic radiation; v„ind is and diffusive escape is 
slow in these models so that the few protons are efficiently 
converted into 7-rays. 

By contrast, in high B models (200-400 ^G), the GeV- 
TeV emission is largely pionic. We find 77 « 0.07-0.14, 
in line with previous use of 77 w 0.1. The electron accel- 
eration efficiency is relatively small, shrinking from 0.02 at 
200 /iG to 0.006 at 300 ^G. The best-fit wind speeds are 150 

- 300 km s"', giving comparable advective and pionic life- 
times for_protonsiiiM82, as expected from the 7-ray bright- 
ness (!Lacki et al.!!201 lb . Diffusion is also relatively quick in 
the B = 200-300 fiG models, probably to compensate for the 
best-fit p = 2.l injection slope and steepen the GeV- TeV spec- 
tra. The B = 400 /iG model, in contrast, has p = 2.3 and no 
diffusion. 

The B = 300 /iG model has the highest S of the best-fit mod- 
els for each B, equal to 49. This is the closest to the Milky 
Way value of 50 - 100, so we adopt it as our fiducial model 
when needed. 

The GeV- TeV results are not significantly affected if we as- 
sume that half of the primary are positrons, instead of all 
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Fig. 5. — The predicted 7-rav luminosity spectram of M82, compared to X -ray (total, stars: 'C appi et"S]|1999l : IMivawaki et al.ll2009t ICusumano et alj|2010l : 
diffuse, square: [Strickland & H eckman 2007) and 7-ray jAbdo et ai]|2010fl : |Acciari et al. 2009) observations. The different components to the emission are 
synchrotron (dash/dotted), positron annihilation (long-dashed/short-dashed), nucleosynthetic 7-ray Hnes (grey lines), bremsstrahlung (short-dashed), IC (dotted), 
and pionic emission (long-dashed). The total emission are the solid lines, black including 77 absorption on the sightline to Earth and grey without it. The plotted 
models show the case for when primary are purely electrons, with no primary positrons. Emission from CCSNe is not included. The X-ray emission mostly 
comes from X-ray binaries and thermal emission from gas and is not included in our fitting. As in Figure|2] note the different energy binnings for the lines. 



of them being electrons. 

Results for MeV j-rays - The main source of MeV 7-ray 
emission is that from Inverse Compton and bremsstrahlung. 
In the fiducial model, the MeV emission is ^ 10% of the 
GeV luminosity (Figure |5] right panel). Because of the small 
scale heights preferred, meaning high density relative to radi- 
ation field, bremsstrahlung is enhanced at expense of Inverse 
Compton. This concentrates luminosity in the bremsstrahlung 
bump at 100 MeV and reduces the power available at the 
lower end of the MeV range. We find that while the 0.1-10 
GeV luminosity of M82 is 1.4 x 10'*'^ ergs s~' in our fiducial 
model, the 1-100 MeV luminosity is only 2. 1 x lO-'^ ergs s"' . 
Surprisingly, we find that synchrotron is the dominant source 
of nonthermal diffuse emission up to 100 keV in our fiducial 
model. This is because the primary electron spectrum extends 
to 7 = 10^ in that model. Considering how little is known 
about PeV CR e^, the actual strength of the synchrotron com- 
ponent is highly uncertain. 

The low B models have stronger leptonic emission, in some 
models erasing the pion bump entirely. The larger scale height 
reduces the importance of bremsstrahlung with respect to IC. 
Even in these models, the bremsstrahlung bump still causes 
the peak of the emission to be at 300 MeV (Figure |5]l. How- 
ever, the relative dominance of IC emission means that uLi, 
is only ^ 3 times smaller at 1 MeV than at 300 MeV in the 
B = 50 /iG model. In the best-fit B = 50 fiG model with no 
primary e"^, the 0.1-10 GeV luminosity is 3.0 x lO'*" ergs s"' 
while the 1 - 100 MeV luminosity is 2.3 x lO'*" ergs s"'. 

The positron annihilation radiation is a possible contrib- 
utor to the 7-ray flux in the 100 - 500 keV range (Fig- 
ure |6]l. In models where the primary are entirely elec- 
trons, the nucleosynthetic positrons are the dominant source 
of positrons. However, even in these models, the relatively 
calorimetric environment means that there is a large rate of 



hadronic positrons cooling all the way down to thermal en- 
ergies: roughly 1/3 of the positron annihilation radiation in 
our fiducial model is from pionic e"^. Thus starbursts have 
enhanced positron annihilation radiation associated with star- 
formation compared to normal galaxies (although non-star- 
formation processes are relatively weaker). In the fiducial 
model, the total positron annihilation energy luminosity is 
3.7 X lO-'^ ergs s"' from positronium continuum and 511 keV 
line emission, plus 7.3 x 10^^ ergs s~' from in-flight annihi- 
lation radiation. Hence the positron annihilation emission is 
only ~ 5% of the total MeV band emission; only near 511 
keV does it actually dominate. 

In models where half of the primary are positrons, 
there is a vast positron annihilation rate, especially in low 
B models which have high ^. These models have com- 
bined positronium and 511 keV line emission luminosities of 
3 X 10^*^ ergs s"' and 5.6 x 10^** ergs s"' for B = 300 fiG and 
50 /iG, respectively, with an additional 2.0 x 10^^ ergs s~' 
and 2.5 x 10^*^ ergs s"' from in-flight radiation. The extraor- 
dinarily intense 511 keV line emission in these models may 
be constrained with future MeV instruments (see Table |2}. 

Finally, the ^^Al and other nuclear lines are a significant 
contributor to the ~ 1 MeV luminosity of M82. The 1.809 
MeV line is the strongest 7-ray line except possibly the 511 
keV annihilation line. However, the photon flux is actu- 
ally small, only ^ 2 x lO'^^cm"^ s"'. As for the nuclear de- 
excitation lines, they are completely negligible in M82 and 
far beyond any proposed detection capability (Table|2). 

High z - The CMB has much less of an effect on a starburst 
Uke M82 than in the Milky Way because the CMB is sub- 
dominant with respect to the high magnetic fields and starlight 
energy density (synchrotron and IC are the dominant loss pro- 
cesses at high electron energy). We find that even at z = 4.9, 
the MeV emission is enhanced only by at most 24% in our 
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TABLE 1 

Best-Fit Models of M82 and NGC 253 



B Prim Best-fit parameters 

e^-' h p fd,ff(3GeV) 7^' v„,„d r, Sff(l GHzf & 

pc Myr km s ' Jy 



M82 



50 


e 


100 


2.1 


00 




1000 


0.025 


0.11 


0.91 


0.44 




e± 


100 


2.1 


00 


10* 


1000 


0.025 


0.11 


0.91 


0.44 


100 


e" 


30 


2.1 


00 


10* 





0.0088 


0.051 


0.41 


0.34 




e± 


30 


2.1 


00 


10* 





0.0088 


0.051 


0.41 


0.34 


150 


e" 


30 


2.1 


10 


10' 





0.025 


0.026 


0.37 


1.9 




e± 


30 


2.1 


10 


10* 





0.025 


0.025 


0.37 


2.0 


200 


e" 


30 


2.1 


1 


lO' 


150 


0.10 


0.017 


0.42 


12 




e± 


30 


2.1 


1 


10« 


150 


0.10 


0.016 


0.41 


12 


250 


e" 


30 


2.1 


1 


10^ 


300 


0.14 


0.011 


0.34 


25 




e± 


30 


2.1 


1 


10* 


300 


0.14 


0.0093 


0.34 


30 


300 


e" 


30 


2.1 


1 


109 


300 


0.14 


0.0058 


0.46 


49 




e± 


30 


2.1 


1 


10* 


300 


0.14 


0.0048 


0.45 


58 


400 


e" 


10 


2.3 


00 


10* 


150 


0.071 


0.023 


0.34 


24 




e± 


10 


2.3 


100 


10' 


150 


0.071 


0.024 


0.34 


23 



NGC 253 Core 



50 


e 


50 


2.1 


10 


10* 


1000 


0.20 


0.19 


0.34 


2.1 




e± 


50 


2.1 


10 


10* 


1000 


0.20 


0.19 


0.34 


2.1 


100 


e" 


15 


2.1 


100 


10* 


150 


0.050 


0.084 


0.20 


1.2 




e± 


15 


2.1 


00 


10* 


150 


0.050 


0.085 


0.19 


1.2 


150 


e" 


15 


2.1 


10 


10* 


150 


0.071 


0.043 


0.18 


3.3 




e± 


15 


2.1 


10 


10* 


150 


0.071 


0.043 


0.18 


3.3 


200 


e" 


5 


2.1 


0.1 


108 





0.10 


0.077 


0.082 


2.6 




e± 


5 


2.1 


0.1 


10' 





0.10 


0.081 


0.076 


2.4 


250 


e" 


5 


2.1 


0.1 


10' 





0.14 


0.050 


0.13 


5.6 




e± 


5 


2.1 


0.1 


10' 





0.14 


0.051 


0.13 


5.5 


300 


e" 


5 


2.1 


1 


10* 





0.071 


0.026 


0.069 


5.3 




e± 


5 


2.1 


1 


10* 





0.071 


0.027 


0.062 


5.3 


400 


e" 


5 


2.1 


1 


10* 





0.10 


0.0071 


0.16 


28 




e± 


5 


2.1 


1 


10* 





0.10 


0.0072 


0.15 


28 



Denotes whether the primaiy electrons and positrons in the model are all electrons {e ) or half electrons and half 
positrons (e^). 
6 I 



Includes both primary electrons and primaiy positrons. 
^ Flux of free-free emission at I GHz observed at Earth. 

fiducial model. In the Z? = 50 jiG model, the MeV enhance- 
ment at z = 4.9 is 65% at 100 keV and less at high energies. 

3.4. NGC 253 Core 

Introduction - NGC 253 (D ss 3.5 Mpc; ' Rekola et all 
|2005) is the only other starburst galaxy detected in both 
GeV and TeV 7-rays, and is s econd only to M82 in in- 
frared brightness for starbursts (ISanders et al.l l2003l l. The 
radio continuum, molecular line, and infrared emission is 
concentrated in a region with a radius between 125 and 
255 pc in radius (e.g.. Turner & Ho 1983; Peng et al. 1996; 
Ulvestad & Antonuccill 19971; lUlvestadll200a ISakamoto et all 



201 lb ; we therefore take the radius of the starburst core to 
be 150 pc. Since the vertical extent of the radio disk im- 
phes a maximum midplane-to-edge scale height of 70 pc 
dUlves tad & Antonucci 1997), we consider scale heights of 
50, 15, and 5 pc. A variety o f measurements give gas 
masses of (2- 5) x 1 0^ M(d (e.g., iMauersber ger et al .lll996l 



Harrison et aO Il999t [Bradford et all l2003l SakamotoetiD 
201 lb . Using a radius of 150 pc gives us a surface density 
of 0.07 — 0.15 g cm~^; we adopt 0.1 g cm"^ here. 

While the total infrared luminosity of NGC 253 as a whole 
is 3.5 X 10'" Lf;^, only half of that comes from the starburst 
itself dMelo et al.ll2002l) . so we adopt a luminosity of 1.7 x 



10'° Lq. Using the lennicua ([1998^ conversion factor and 
translating this to a SalA IMF, we get a star-formation rate of 
2.3 M0 yr"' and a supernova rate of 0.020 yr"' (compared 
to the c anonical 0.03-0.1 yr"'; again, see the summary in 
iLacki e t al. 2011). 

Like M82, almost nothing is known about NGC 253's MeV 
emission from observations. It is not even detected in 10 - 
50 keV X-rays by Swift-BAT (ICusumano et all 1201 Oh . In- 
terestingly, OSSE did report a detection for NGC 253 in the 
100 keV range, which iGoldsh midt & Re phaehl d 19951) inter- 
preted as IC emission from the halo of NGC 253. How- 
ever, this detection implies a ^ MeV luminosity many times 
greater than the GeV luminosity detected by Fermi. Fur- 
thermore, NGC 253 was not detected in this energy band 
with the more sensitiv e Swift-BAT, so we consider it spurious 
(Cusumano et al.ll20T 0). although it could have been from a 
transient source. It is interesting, though, that if the detection 
were typical of starburst galaxies, then starbursts would make 
up the MeV background. OSSE have provided very weak 
const raints on the continuum and line emission from NGC 
253 (Bhatt acharya et al.l ll994'). In the two-year Ferm i-LAI 
catalo g, NGC 253 is detected down to 300 MeV ( Nolan etalJ 
I2OI2I) . A number of theoretical models o f the nonther- 
mal emission predict the MeV spectrum (e.g.. lPaglione et al] 
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TABLE 2 
Gamma-ray Line Strengths 



Line^ 




Source'' 




M82 






NGC 253 Core 








LowB*^ 


Hij 


jhS'' 


LowB'= 


Hig 


hB' 








Luminosity 


Photon Flux 


Luminosity 


Photon Flux 


Luminosity 


Photon Flux 


Luminosity 


Photon Flux 




MeV 




ergs s"' 


ph cm"^ s"' 


ergs s"' 


1 —9-1 

ph cm s 


ergs s"' 


ph cm"- s"' 


ergs s"' 


1 -7-1 

ph cm s 


annih. 


0.511 


CR (hadr.) 




2 X lO" 




1 X 10"'" 


3 X lO^'' 


3 X 10"' 




2 X 




1 X 10"'" 


3 X 103'^ 


2 X 10"' 






CR (prim. £■+) 


3 X 10^' 


3 X 10"'' 


1 X 103* 


1 X 10"'' 


2 X 103' 


1 X 10"'' 


6 X 103^ 


5 X 10"* 






Nucleo. 


2 X 10^^ 


1 X 10"^ 


2 X 10" 


1 X 10"* 


5 X 103'' 


4 X 10 ^ 


5 X 103* 


4 X 10"' 


^''Al decay!^ 


1.809 


Nucleo. 


7 X lO" 


2 X 10"^ 


7 X 10" 


2 X 10"* 


2 X 103' 


5 X 10"' 


2 X 103^ 


5 X 10"' 


^^Fe decay^ 


1.173 


Nucleo. 


5 X 10^'' 


2 X 10"' 


5 X 10^6 


2 X 10"' 


2 X 103'' 


6 X 10"'" 


2 X 103^ 


6 X 10"'° 




1.333 


Nucleo. 


5 X 10^'' 


2 X 10"' 


5 X 10^' 


2 X 10"' 


2 X 103'' 


6 X 10"'" 


2 X 103*^ 


6 X 10"'° 


^Ti decay^ 


1.157 


Nucleo. 


1 X 10^'' 


4 X 10"' 


1 X 10" 


4 X 10"' 


3 X 103* 


1 X 10"' 


3 X 103'5 


1 X 10"' 


^He de-ex. 


0.429 


CR (broad) 


7 X 10^' 


7 X 10"'^ 


4 X lO^'^ 


4 X 10"'3 


2 X 1032 


2 X 10"'3 


1 X 103^^ 


1 X 10"'3 


^He de-ex. 


0.478 


CR (broad) 


9 X 10^' 


8 X 10"'" 


6 X 10" 


5 X 10"'3 


2 X 1032 


2 X 10"'3 


1 X 1032 


1 X 10"'3 


'^C de-ex. 


4.438 


CR (narrow) 


7 X 10" 


6 X 10"'3 


6 X lO^" 


5 X 10"'2 


1 X 103" 


1 X 10"'^ 


2 X 103" 


2 X 10"'2 






CR (broad) 


2 X 10'-* 


1 X 10"'2 


1 X 10" 


9 X 10"'2 


3 X 103" 


3 X 10"'2 


2 X 103" 


2 X 10"'2 


de-ex. 


6.129 


CR (narrow) 


4 X 10" 


2 X 10"'3 


3 X 10'" 


2 X 10"'2 


7 X 1033 


5 X 10"'3 


1 X 103" 


8 X 10"'3 






CR (broad) 


6 X 10" 


4 X 10"'3 


4 X 103" 


3 X 10"'2 


1 X 103" 


1 X 10"'2 


9 X 10" 


6 X 10"'3 


'*'0 de-ex. 


6.916 


CR (nan'ow) 


1 X 10" 


6 X lO""* 


7 X 1033 


4 X 10"'3 


2 X 1033 


1 X 10"'3 


2 X 10" 


1 X 10"'3 






CR (broad) 


2 X 10" 


1 X 10"'3 


1 X 103" 


6x 10"'3 


4 X 1033 


2 X 10"'3 


2 X 10" 


1 X 10"'3 


"=0 de-ex. 


7.115 


CR (nan'ow) 


1 X 10" 


6 X 10"'" 


8 X 1033 


4 X 10"'3 


2 X 1033 


1 X 10"'3 


2 X 10" 


1 X 10"'3 






CR (broad) 


2 X 10" 


1 X 10"" 


1 X 103" 


6 X 10"'3 


4 X 1033 


2 X 10"'3 


2 X 10" 


1 X 10"'3 


'^C de-ex. 


15.10 


CR (narrow) 


1 X 10" 


4 X 10"'5 


1 X 1033 


3 X 10"'" 


3 X 1032 


8 X 10"'^ 


4 X 1032 


1 X 10"'" 






CR (broad) 


5 X 10" 


1 X 10"'" 


3 X 1033 


8 X 10"'" 


1 X 1033 


3 X 10"'" 


7 X 1032 


2 X 10"'" 



^ "Annih." abbreviates annihilation and "de-ex." abbreviates de-excitation. 

'' "Nucleo." stands for 7-ray lines produced by short-lived radioisotopes produced directly by stellar isotopes, while "CR" stands for 7-ray lines powered by CRs interacting with the ISM. 

For de-excitation lines, "narrow" stands for the line component produced by the ISM atoms after being hit by a CR, while "broad" stands for the Doppler-shifted line component emitted by 

the CRs after hitting an ISM atom. For the positron annihilation line, "hadr." stands for hadronic positrons and "prim stands for primary positrons. 

^ Values for the B = 50 /iG model of M82 with no primary positrons, except for the "CR (prim, e^)" row, which refers to the model with primary electrons and positrons. 

'' Values for the B = 300 i^lG model of M82 with no primary positrons, except for the "CR (prim, e^}" row, which refers to the model with primary electrons and positrons. 

^ Values for the S = 50 i^lG model of NGC 253 core with no primary positrons, except for the "CR (prim, e^)" row, which refers to the model with primary electrons and positrons. 

^ Values for the B = 400 fiG model of NGC 253 core with no primary positrons, except for the "CR (prim, e^)" row. which refers to the model with primary electrons and positrons. 

^ The nucleosynthetic 7~ray lines have exactly the same luminosities and fluxes in the high B and low B models. 



1996 : iDomingo-Santamarfa & Torresl 120051: iRephaeli et al.l 



2009 ) 



Constraints from radio and GeV-TeV j-rays — We find gen- 
erally similar results for NGC 253's starburst core^ as we do 
for M82, which is not surprising given the similar physical 
conditions in these starbursts (see Table[T}. 

The radio spectrum for the core of NGC 253 is even flatter 
than that of M82 (see Fig.|7|. In addition to forcing p to 2.1, 
it forces the best-fit scale heights to be even smaller than for 
M82: 15 pc for B = 100- 150 pc and 5 pc for higher B. This 
is because of the flattening effect of bremsstrahlung a nd ion- 
izatio n losses on the electron spectrum (c.f., Thompso n et alj 
|2006a). Of course, scale heights of only 5 pc for NGC 253's 
starburst are not physically realistic. Another possibility, here 
unexplored, is that CRs are preferentially interacting with gas 
of higher density than the mean: for example, if they only 
interact with the dense cores of molecular clouds. If true, 
this would be the opposite of the apparent situation of the 
Galactic Center, where the low 7-ray luminosit y suggests that 
CRs av oid the dense cores of molecular clouds (ICrocker et al.l 
1201 Ibl) . 

As with M82, the GeV emission in low B models is mainly 
leptonic and in high B models it is mainly hadronic. The TeV 
emission is pionic, requiring the primary spectrum to cut 
off at 7m"'x in low B models. At B = 50 /iG, the fit is poorer 
than in the other models, suggesting that Z? > 50 /iG in NGC 
253's core as well (x^ = 110 for 50 fiG instead of 70-80 
for higher B), although the difference in fit quality is not as 

* Throughout the paper, we use the tenns "core" and "starburst" for NGC 
253 interchangeably, to refer to the modelled regions within 150 pc of its 
center 



great as in M82. We find that ry approaches ^0.1 for B > 
200 /iG, as expected. However, we find small wind speeds of 
150 km s~' in 100- 150 /iG models and no winds in models 
with higher B. This may be because h is so small in the high 
B models, so that even a slow wind rapidly advects CR nuclei 
out of the starburst; if instead both /; and («) are larger in 
reality, more rapid winds would be allowed. On the other 
hand, rapid diffusion (0.1-1 Myr at 3 GeV) is required in 
the best-fit models with B > 200 /iG, probably to steepen the 
pionic 7-ray spectrum from F « 2.1 to the observed 2.3. 

None of the best-fit models have Milky Way-like values of 
S. In fact, all of them with B < 300 /iG have (5 « 1 - 6, which 
would indicate significant power in leptons. The B = 400 /iG 
model has 6 = 28, about a third of what we would expect in 
the Milky Way; we adopt it as our fiducial model. 

Like M82, the GeV- TeV results are not significantly af- 
fected by assuming a mixture of primary electrons and 
positrons instead of primary electrons. 

Results for MeV ^-rays - Much like M82 models, the pi- 
onic and bremsstrahlung emission concentrates power at 100 
MeV and above, and leptonic emission at lower energies is 
relatively weak, especially if B is high. The small 5 even in 
our fiducial model increases the amount of leptonic emission. 
However, the small scale heights (and thus, high gas densities) 
selected pushes that leptonic emission into bremsstrahlung. 
Thus, the 1-100 MeV emission is even weaker relative to 
the pionic emission than in M82. In our fiducial (B = 400 /iG) 
model, we find a 1 - 100 MeV luminosity of 2.1 x lO-'^ ergs s"' 
and a 0. 1 - 10 GeV luminosity of 1 .4 x lO"**' ergs s"' . The vLjj 
luminosity at 1 MeV is just 2% that at 1 GeV. 

The low B models of course have stronger MeV emission 
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Fig. 6. — The continuum positron annihilation radiation (thin solid), both 
positronium continuum (black) and in-flight annihilation (grey), for the best- 
fit B = 300 /xG model of M82 where half of the primary are positrons. The 
contributions from hadronic e+ are short-dashed; nucleosynthetic e+ from su- 
pemovae are long-dashed, and the possible primary component is dotted. The 
supposed primary positron component dominates the low energy e+ popula- 
tion, although hadronic in-flight radiation dominates over 30 MeV. Nucle- 
osynthetic e+ dominate over the hadronic e+ below 1 MeV. For comparison, 
the total emission from M82 in this model is shown as the thick solid light 
grey line. Also shown in the total emission for the fiducial B = 300 fj,G model 
with no primary positrons (thick dot-dashed grey line); the low energy con- 
tinuum shape is different because a higher 7n™ is preferred slightly by our 
fitting, resulting in a much higher synchrotron contribution. 



NGC 253 Core 
B = 400 /J.G 




100 



1/ [GHz] 



Fig. 7. — Predicted radio spectrum of NGC 253's starburst in the best-fit 
B = 400 fiG model with no primary positrons. Plotted Hues are total radio 
emission (solid), total synchrotron emission (long-dashed), synchrotron from 
secondaries (short-dashed), synchrotron from primaries (dotted), and free- 
free emission (dash-dotted). Plotted data are the interferometric data com- 
piled in Williams & Bower (2010) (open circles). The bottom panel shows 
the residuals to the radio data. 

redshifts. At z = 4.9, the 100 keV emission of the fiducial 
model is enhanced by 60% and the 1 MeV emission in en- 
hanced by a third. The CMB has a greater effect in the 50 fiG 
model, enhancing the 100 keV emission by 90% at z = 4.9 and 
the 1 MeV emission by 50%. 



(Figure|8]l: with no primary positrons, the 1 - 100 MeV lumi- 
nosity is 2.3 X lO''" ergs s"' and the 0.1-1 GeV luminosity is 
3.0 X 10"*" ergs s"'. 

Like M82, the higher proton calorimetric efficiency en- 
hances the positron annihilation rate per unit star-formation, 
with about 1 /3 as many hadronic positrons as nucleosyn- 
thetic positrons in our fiducial model. In the fiducial model, 
positrons annihilating near rest contribute 1 .4 x 10^^ ergs s"' 
through positronium continuum and 511 keV line emis- 
sion, and relativistic positrons contribute an additional 5.5 x 
10^^ ergs s~'. Thus, positron annihilation does not signif- 
icantly increase the total MeV luminosity, although it may 
dominate near 511 keV. Models with primary positrons have 
much higher annihilation luminosities: for B = 400 /iG, an- 
nihilation near rest produces 1 .2 x 10^"* ergs s"' and in-flight 
annihilation produces 1.1 x 10""^ ergs s"' ; for B = 50 /iG, anni- 
hilation near rest produces 2.7 x 10-''^ ergs s"', while in-flight 
annihilation produces 1.2 x 10^' ergs s"'. 

We list the strengths of the 7-ray lines in Table |2] The 
strongest is the ^^Al 1.809 MeV line, with a luminosity of 
2.7 X 10^^ ergs s"'. None of the lines are expected to con- 
tribute significantly to the MeV luminosity of NGC 253, or be 
detectable with planned instruments. Note that there is signif- 
icant star-formation outside the starburst of NGC 253 and in 
the main disk, w hich will enhanc e the unresolved line fluxes 
by a factor < 2 (Melo et al. 2002). 

High z - As with M82, the high magnetic field and starlight 
energy densities limit the influence of the CMB to very high 



3.5. Brief Summary 

The starbursts of M82 and NGC 253 are more luminous 
per unit star-formation than the Milky Way at GeV ener- 
gies, but this is due to strong pion losses for protons, lead- 
ing to a prominent pion bump in the "high B" models. At 
present, we have essentially no MeV data for starbursts, so 
we cannot empirically rule out leptonic "low Z?" models with 
strong MeV emission, though these require strange parameter 
choices. Compared to the high-B starburst models, the Milky 
Way has a higher MeV/GeV ratio, due to relatively efficient 
conversion of primary e~ energy into IC and bremsstrahlung 
radiation. If it is placed at high z, the Milky Way's MeV- 
GeV spectrum becomes flat from increased Inverse Compton 
losses off the CMB. We therefore a priori do not expect ei- 
ther starbursts (with their prominent pion bumps) or normal 
galaxies (which are relatively dim at GeV to begin with) to be 
particularly strong MeV sources. 

4. MeV EMISSION FROM INTERGALACTIC CASCADES 

Even if the 7-ray spectra of star-forming galaxies were 
purely pionic, MeV emission is produced during the propa- 
gation of 7-rays from distant star-forming galaxies to Earth 
through the cascade process. About half of the star-formation 
in the Universe occurs before z « 1 . At these distances, the 
Universe is opaque above 100 GeV from pair production pro- 
cesses (7+7 —> e"'" + e~) on the in frared Extragalactic Back- 
ground Light (EB L; IStecker et al1[T992' Gilmore et al. 20091 
iFinke et al]|2010l) . The have a typical energy comparable 
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Fig. 8. — T he predicted 7-ray luminosity spectram of NGC 253's starburst, compared to X-ray fCappi et al.l ll999l : IBauer et al.ll2008 l) and 7-ray <Abdo et aP 
l20101tlAcero et al. 2009) observations. The different components to the emission ai'e synchrotron (dasfi/dotted), positron annihilation (long-dashed/short-dashed), 
nucleosynthetic 7-ray lines (grey lines), bremsstrahlung (short-dashed), IC (dotted), and pionic emission (long-dashed). The total emission are the solid hnes, 
black including 77 absorption on the sightline to Earth and grey without it. The plotted models show the case for when primary are purely electrons, with 
no primary positrons. Emission from CCSNe is not included. The X-ray emission mostly comes from X-ray binaries and thermal emission from gas and is not 
included in our fitting. As in Figure|2] note the different energy binnings for the different lines. 



to the original 7-ray, and cool mainly by Inverse Compton 
scattering of the CMB. 

The typical rest-frame energy of upscattered CMB light is 
then 



;32 MeV(l-i-z) 



100 GeV 



(8) 



where z is the redshift where the scattering takes place. Thus 
intergalactic cascades shift energy from TeV bands, where pi- 
onic emission can be important, to MeV bands. Normal star- 
forming galaxies have a pionic spectrum that goes as E"^-^, 
with little 7-ray emission above 100 GeV, so they will con- 
tribute little cascade emission. Starburst galaxies, on the other 
hand, are observed to have hard E"^^ spectra between GeV 
and TeV energies. Nearly 1 /3 of their > GeV luminosity is 
above 100 GeV and can be affected by the cascade process. 

Recentlv. lBroderick et a l. (2012) have proposed that plasma 
instabilities alter the evolution of cascades from bright TeV 
sources like blazars. These instabilities would effectively stop 
before they could radiate IC, and quench all of the cas- 
cade emission. However, the TeV luminosities required for 
this process are much higher than from M82 and NGC 253. 
IChang et alJ (120^1 2) conclude that only HyperLIRGs would be 
affected by the plasma instability by assuming that the IR/7- 
ray ratio is the same as in M82 and NGC 253. We note that 
purely proton calorimetric galaxies would be expected to have 
IR/7-ray ratios 2-3 times higher than M82 and NGC 253's 
starburst, and so TeV emission from brighter ULIRGs (which 
are numerous at z w 2) may not cascade. However, since most 
star-forming galaxies are relatively faint, we ignore the effects 
of plasma instabilities. 

4. 1 . Calculation of Cascade Radiation 

Suppose the star-forming galaxy lies at redshift Zs and emits 
a spectrum dN^'''^^ /dE' over some time interval. Then the 



cascade spectrum observed at redshift Zobs can be calculated 

as 

(E^ ,E ) 



^^„-(4,Zobs)= / dz 



7 



/o dzdE^' dEi;- " 

X fiife exp(-T-y-y(£'^ , Zob.s , z))dE'^ (9) 
where = (\+z)E^ is the cascade-frame energy, E'^ = 
(1 +Zohs)E'y is the source-frame energy, and E^ is the Earth- 
frame energy. The key ingredients in this calculation are 
dNe/{dE'Jdz), the spectrum of pair generated in the red- 
shift step z to z- dz from both the primary 7-ray spec- 
trum and the cascade 7-ray s at previous redshift steps; 
dQic/dEj(Ee,E^), the IC emission at Ey from one CR of 
energy E/, fic, the IC cooling time for e^; and Tyj(Ej,Zohs,z), 
the 77 optical depth to a photon observed at E-^ from redshift 

Zobs to Z. 

Accurate source functions for pair production are given in 
Aharonian et al. (1983) and Bottcher & Schlickeiser (1997). 
We consider the emission from star-forming galaxies at red- 
shift Zs separately. We then calculate the cascade emission in 
redshift steps, each from z, to z, - dz- At each step, the input 
photon spectrum (in the cascade rest frame) is 

dN,.. 



dNy _ ^/A?;™"^ 1+z. 



dEL'^ 



dE'y 



1+z 



exp(-T.^^(Zobs,z))+' 



dE'^ 



■(E''z) (10) 



which includes both the attenuated primary spectrum and 
the attenuated cascade spectrum from previous redshift steps. 
Within the redshift step, the UV-IR background effectively 
forms a screen to the 7-rays with optical depth: 



A'. 



77 ■ 



:T..^^(£'",z/-Az,-,Zi) 



So the generated pairs are calculated as 

dN, f dNy dQe(E'^) dt 1 - exphAr^^,] 



dzdE'- 



dE^ 



dE'J dz 



At, 



(11) 



dE'' (12) 
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where dQ(E!^)/dE'J is the rate at which pair at energy 
are produced by a photon of energy E^, dt/dz puts the red- 
shift step in terms of time. The exponential term guarantees 
energy conservation even if the Universe is opaque through 
the redshift step; in the limit of large r the dz{dt/dz){\- 
exp(-r)) /(r) expression reduces simply to the timescale for 
photon- photon annihilation.^ [ Bottcher & Schlickeiserl (1 1997b 
give the lAharonian et al expression for the rate of 

pair production per high energy photon of energy einieC^ as: 



dQAE'') 3 car 



dE'j 



16 m-jC^ 



d£2- 



1 finback 



2(2£i£2- 1)£ 



(£2)^ d£2 

\-i> 

( 4£27(£i-7) 



8£i£2 



7(£l-7) 



1 

£\£2 



T(£\-l) 



a 3) 



where £2 is the cascade-frame energy of a background photon 
in units of m^c^ and dn\,.^ck/de2 is the cascade-frame number 
spectrum of background photons. 

The rate of Inverse Compton cascade emission from each 
pair is calculated as 



dQ 



ic 



dEi: 



■■cE' 



' ,dn(e',z) I, I 
dt —jp—cric(E^,e ,E,) 



(14) 



where dn(e' ,z)/de' is the rest-frame backgrou nd radiation 
spect rum at redshift z and rest-frame energy e' (ISchlickeiserl 
I2OO2I) . The IC cross section is equal to 



ciciE^je ,Ee)-- 



3aT 



2^1n^+(l + 2^)(l-^) + 



(15) 

with q = E y/[r,(E,- E^] (1 /(47^) < ^ < 1) and = 
4e'7/(mc2) dSchlickeisedl2QQ2h . 

In the intergalactic medium, electrons and positrons are 
mainly cooled by Inverse Compton losses and the adiabatic 
losses from the expansion of the Universe (again ignoring any 
plasma instability losses). The electron lifetime to Inverse 
Compton emission is 



tic(E'J) = E', 



'E''^dE'' 
^ dEli ^ 



(16) 



while the loss time due to the expansion of the Universe are 
simply 

tH = [Ha^/^Mil + zf + nAT'- (17) 
The final lifetime is then just found as ty^^^ = fj^' +f^'. 

4.2. Results of Cascade Calculation 

The effects of cascades on the spectra of Milky Way-like 
galaxies is in fact negligible, as expected. However, cascades 
do have a small but noticeable effect on the spectrum of our 
high-B M82 spectra, as shown in Figure |9] Cascade emis- 
sion makes up ~ 10% of the 100 MeV - 10 GeV emission of 
starburst galaxies at z = 2. There is a spectral downturn be- 
low 100 MeV, however the intrinsic emission of these models 

' The term comes from the time-averaged mean of the photon flux within 
the redshift step, which goes as ^ e"^ dr / ^ dr. 



Note that formulas 26 and 32 of lBiittcher & Schlidceiserl 09971) are for 
electrons only, not for both electrons and positrons (Bottcher, personal com- 
munication; Aharonian & Khangulyan, personal communication). 




1 2 3 4 5 
log,„ [E; (MeV)] 

Fig. 9. — The effects of intergalactic cascades on M82's apparent rest- 
frame spectrum shape, for the fiducial model placed at z = 2. The grey dotted 
line is the unpropagated input spectrum (after absorption within M82 itselt7 
the grey solid hne is the cascade emission, and the solid black line is the 
attenuated plus cascade emission. For comparison, the unpropagated pionic 
spectrum is shown as a dashed grey line. The cascade spectrum sticks out of 
the pion bump, but is less strong than the leptonic emission. 

falls even faster because of the pion bump. Thus, the fraction 
of emission due to cascades peaks near 10 MeV, where it 
makes up about a third of the observed emission from star- 
bursts. 

We conclude that cascade emission alters the MeV spec- 
trum of starburst galaxies at high z, but nowhere near enough 
to fill in the pion bump by itself. 

5. CONSTRAINTS ON COSMIC MeV EMISSION FROM 
STAR-FORMING GALAXIES 

5.1. Using the templates to predict the background 

The number density of 7-rays observed at energy E^ is 



dn 
dE, 



dt- 



d^qo 



(18) 



where dt = t/z/[-Ho(l + z)\/^A + ^Mi)-+z)^} describes the 
lookback time, Vcom is the comoving volume, and 
d^qcom/dVcomdE^ is the injection rate of 7-rays per unit vol- 
ume per unit observed energy. The number density is con- 
verted into a photon flux by multiplying by c/(4-tt). This gives 
us 

c d qQQYn 



d^ 
dE. 



-(£7) = 



dz- 



"ATiHail+zWMil+zy + ^KVl^ dVcomdE^' 

(19) 

This expression can also be derived using the photon number 
luminosity distance. The energy flux per (natural) log bin of 
energy is i^/^ = E^d^/dEj. 
We calculate the injection rate of 7-rays as 



SFR, 



model 



dEi 



using p^™ is the star-formation rate per unit volume, f(z) is 
the fraction of star-formation from galaxies like the model 
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galaxy, dQ j dE'^ is the total source-frame photon generation 
rate in the model galaxy (after cascades), and SFRmodei is the 
star-formation rate of the model ga laxy. The co moving star- 
formation rate comes from Hopkins "& BeacomI (j2b06). 

To compare the predicted background to the observed back- 
ground, we link together fits to the X-r ay background and the 
GeV background. At X-ray energies, lAiello et all (l2008h fit 
the X-ray background with a double power law, while at GeV 
energies, lAbdo et aP (1201 Oah fit the GeV background with a 
single power law. We connect the two with a power law (lin- 
ear fit in log E and log vl^) at intermediate energies, which 
we take to be 3 to 20 MeV by visual inspection. We use the 
resulting function to represent the strength of the X-ray-7-ray 
background: 



vL, 



Q.lQlEl 



Eke 



1.32 



1.58 



1.45 



29.99 

^keV 



£'keV 

29.99 



0.303 



2.88 



3000 

£'keV 

100000 



(£ < 3 MeV) 



(3 MeV < £ < 20 MeV) 



-0.41 



(E > 20 MeV) 



in units of keV cm s 'sr for a photon of energy E. 



(21) 
with 



£kev = £/keV (AjelloeLalJ|2008; Abdo et al. 2010a). The 
function is plotted in Figure [T] (and Figures [TO] [TT] \l7\ as 
the solid dark blue line. 

5.2. The question of the starburst fraction of the cosmic SFR 

A large uncertainty is the fraction of cosmic star-formation 
history that occurs in starburst galaxies. Indeed, discussions 
of the starburst contribution to the 7-ray background must 
be clear about what "starburst" even means for these calcu- 
lations. For this paper, we simply define starburst as having a 
similar 7-ray SED (and thus, CR propagation environment) as 
M82 and NGC 253. In particular, it means (1) similar 7-ray 
to SFR ratios as M82 and NGC 253's core and (2) hard GeV - 
TeV 7-ray spectra, from winds and pionic losses. At t h e high 
end of starburst fraction estimates. Loeb & Wax maiil (l2006b 
and Thompson, Ouataert, & Waxman (2007) argued that most 
of the star-formation of the Universe was in proton calorimet- 
ric galaxies, and thus, the 7-ray backgrou nd from starbursts 
is high. Further down in starburst fraction. ISteckerl (l2007b ar- 
gued that most of the star-formation at z < 1 was in normal 
galaxies, implying a small 7-ray background from starbursts. 
At the low end, iStecker & VentersI (1201 lb have recently used 
semi-analytic models which argue that the fraction of cos- 
mic star-formation in pure merger-driven starbursts is very 
low, ~ 1-2% at z = 0, 3% at z = 1, and - 5% at z = 4 
(iHopkins et a l.l l20101 

There are potentia l issues wit h each of 

these approaches. iLoeb & WaxmanI (l2006h and 

[Thompson. Ouataert. & WaxmanI (|2007|) m^gue that nor- 
mal galaxies have ~ 10 times more gas at high z than at 
low z; therefore they should be ^ 10 times more efficient at 
pionic losses and nearly proton calorimetric. However, at 
sufficiently high gas surface densities, advective losses may 
be strong enough to negate this effect. In fact, the Galactic 
Center has similar gas surface densities (Sg « 0.03 g cm"^), 
but is in fact fainter in 7 -rays per unit star-fo rmation than 
the surrounding Galaxy (ICrocker et al.l 1201 lab). This is 
apparently because of a strong wind in the region and also 
because CRs do not appear to be fully sampling the gas 



( ICrocker et al.l 1201 Ibl) . Therefore, normal galaxies at high 
z may form stars quickly and have large amounts of gas 
without having M82-like 7-ray spectra (although the pionic 
7-ray spectrum will likely be hard because both pionic and 
advective losses harden the spectrum). However, galaxies 
at high z appear to li e on the far-infrared correlation (e.g., 
lAppleton et al. 2004: iGam & Alexanded 120091: ISargent et all 
I2OIO; Mao et al. 201 1), whereas the Galactic Center does not 
(Crocker et al. 2011a), which is consistent with CRs in high 
z galaxies interacting with average density gas and radiating 
before escaping (Lacki et al. 20 10). 

The predictions of Stecker & VentersI (|201 1|) ar e based on 
the semi-analytic models of .Hopkins et al.l (1201 Oh (HIO). In 
these models, galaxies have a normal star-formation mode 
most of the time, but have a finite probability of merging; 
when they merge, they enter into a brief (^100 Myr) starburst 
mode of star-formation whose evolution is gui ded by gas hy- 
drodyn amics simulations. The simulations of iHopkins et al.l 
(1201 Ol) do not include processes like stellar bars which may 
feed the lower luminosity starbursts. For example, M82 (the 
starburst with the largest 7-ray flux) is interacting with its 
neigh bor M8 1 and has a compact starburst driven by a stellar 
bar (T elesco et 31.1119911) . but is not yet merging, so it might 
not be counted by HIO. Also, a normal galaxy with a high 
enough gas density and weak enough advective losses will be- 
come proton calorimetric and will have a starburst-like 7-ray 
spectrum, even if it is not formally starbursting. Thus defining 
starbursts purely as ^ 100 Myr long mergers may underesti- 
mate the contribution of starburst-like 7-ray emitters. It is 
nonetheless likely that the strongest compact starbursts like 
Arp 220, which arise from true mergers, are a small fraction 
of the cosmic star-formation rate. 

Observationally, a number of lines of evidence point to a 
starburst fraction of ~ 10%, in terms of star-formation mode. 
The total core-collapse supernova rate withi n 10 Mpc is ^ 
0.3-2 yr"' (e.g., Kochanek et al. 2008; Kis tler et al.ll2011l : 
iLi et aLll201 It [Horiuchi et al.ll201 II) . The compact starbursts 
in M82, NGC 253, and NGC 4945 (D w 3.9 Mpc) have a to- 
tal supernova rate o f ^ 0.12 yr"' from their FIR luminosity 
(ISanders et al.l l2003h . Further out, there are nuclear starbursts 
in NGC 463 \{D = 1.1 Mpc) with Tsn = 0.01-0.06 yr"' (Colli 
[19991) . and smaller nuclear starbursts in M51, M83, Maffei 2, 
NGC 2903, NGC 4736, and IC 342 ( Kennicutt.199 8). Conser- 
vatively assuming a total starburst supernova rate of 0.1 yr"' 
within 10 Mpc, we find a starburst fraction of ^ 5-30%. 
Radio and IR luminosity functions can be fit with two pop- 
ulations of star-forming galaxies, identified as normal galax- 
ies and starbursts ( Yu n et al.. 2001.) . These luminosity func- 
tions indicate that starbursts make up ^ 10- 20% of the star- 
forma tion rate at z ~ (e.g., I Yun et al1l2001t iBothwell et alJ 
I2OI Ih . Whether these star bursts are all compact a nd 7-ray 
bright like M82 is unclear. iRodighiero et al.l (1201 1 *) find that 
only 10% of the cosmic star-formation is in starbursts that 
lie off the "main sequence" of galaxies, a rel ation between 
the st ar-formation rate and stellar mass (see also 'Sargen t et alJ 
I2OI2I) . Hopkins & Hemquist (2010) derive a starburst frac- 
tion of 5- 10% by using the surface brightness profiles 
of spheroidal galaxies. In addition to the compact starbursts 
at z ~ 0, there are the submillimeter galaxies (SMGs) ob- 
served at z > 2. There is still some debate whether SMGs 
are simply disk galaxies with very high star-formation rates 
or mergers, but they appear to make up 2 - \0% of the cosmic 
star formation (e.g., Chapman^t al.l 120051: iDve et al.l 120071: 
iMichalowski et al.ll2010hT With vast amounts of gas, SMGs 
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Fig. 10. — The effects of starburst fraction on the predicted 7-ray back- 
ground. Starbursts (fiducial M82) are in black and nornial galaxies (fiducial 
Milky Way) are in grey. The long-dashed line is the low starburst fraction his- 
tory, the solid line is the fiducial medium starburst fraction, the short-dashed 
line is the high starburst fraction, and the dotted line is if all of the star for- 
mation is in starbursts. The nonnal galaxy contribution is only different for 
the high starburst fraction, whereas the starburst contribution can vary by a 
factor ~ 30 between the different starburst fractions. The observed 7-ray 
background is shown in blue. 

are likely to have efficient CR proton losses, although their 
large sizes means that the physical conditions (and therefore 
leptonic emission) could be different than M82. 

In this work, we consider three cases; (1) the low starburst 
fraction case, calculated by integra ting over the Schechter 
functions given in HIO (similar to the lStecker & Ventersl201 ll 
calculation), rising from 1 .6% at z = to 3.8% at z = 2 to 5.1% 
at z = 4 and beyond "; (2) a h igh starburst fraction case as 
given in Thompson, Ouataert. & Waxman (2007); and (3) our 
fiducial medium starburst fraction case, which is 15% at all 
redshifts, since that is roughly the geometric mean of the low 
and high cases and it is fairly close to the ^ 10% that seems 
to hold at low z- However, without further consideration of 
the gas evolution in the normal galaxies, we note the fiducial 
model is not necessarily accurate. 

5.3. The predicted MeV - GeV j-my background 

Using the cosmic star-formation rate, assumptions about 
the contribution of starbursts to the cosmic star-formation his- 
tory, and our Milky Way and M82 templates, we can tenta- 
tively predict the 7-ray background with eqn. [19] We list our 
results for different galaxy spectral templates and starburst 
fractions in Table [3] 

In general, we find that starburst galaxies have a lot of ab- 
solute power per unit star-formation at GeV energies, but nor- 
mal galaxies have much relative power per unit star-formation 
at MeV energies (Figure [TJ. Using the fiducial/medium star- 
burst fraction history, we find that normal galaxies contribute 
~ 6 - 15% of the 0. 1 - 10 GeV background. The normal galax- 

" Since the HIO luminosity functions are not designed for use at z > 4, 
we use the z = 4 starburst fraction at these redshifts. 



ies then contribute ^ 2-4% of the 10 MeV background, < 2% 
of the 1 MeV background, and - 0.1-0.2% of the 100 keV 
background (neglecting discrete hard X-ray sources in these 
galaxies). Note that this drop-off in the 7-ray background 
contribution by two orders of magnitude is not largely due 
to the pion bump: as can be seen in Figure [T] ;//^ for fidu- 
cial normal galaxies only drops off by a factor ^ 2 from 0. 1 
to 100 MeV. The discrepancy arises instead because the 7-ray 
background is steeply falling with energy. 

The starburst contribution to the 7-ray background depends 
on the uncertain contribution of starbursts to the cosmic star- 
formation rate (§ 15.21 and Figure [TOt and, for MeV energies 
in particular, the magnetic field strengths of the starbursts. 
We find the starburst contribution at 1 - 100 GeV is only 

0.7-8% with the low starburst fraction, --4-41% with 
the medium starburst fraction, and from ^ 16-200% with 
the high starburst fraction (Table |3}. However, the starburst 
contribution rapidly falls off at lower energies because of the 
prominent pion bump: with the medium starburst fraction and 
the fiducial M82 spectral template, the contribution falls to 
2.8% at 100 MeV, 0.4% at 10 MeV, 0.1% at 1 MeV, and 
0.01% at 100 keV (again, ignoring discrete sources within the 
galaxies). Note that the pion bump in high B models does 
not mean that starbursts are MeV dim per unit star-formation 
compared to normal galaxies; it arises because starbursts are 
Ge V-TeV bright through their higher proton calorimetry frac- 
tion. Of course, in the leptonic B = 50 //G models, the star- 
burst contribution at MeV energies is much higher: 8 - 13% at 
100 MeV, 4-6% at 10 MeV, 1-2% at 1 MeV, and 0.08-0.3% 
at 100 keV or the medium starburst history. 

The 7-ray lines, whether from CCSNe, e""" annihilation from 
nucleosynthetic positrons, or unstable isotopes in the ISM, 
do not constitute the majority of the MeV contribution from 
star-forming galaxies; continuum processes are the dominant 
source of MeV luminosity of star-forming galaxies. The lines 
do make up ^ 1 /4 of the star-forming galaxy contribution at 
0.5 MeV and > 10% from -0.2-1 MeV. The fraction of the 
MeV emission from star-forming galaxies from the transient 
CCSNe emission peaks at — 17% at 0.4 MeV. The exception 
is in non-standard low B models of starbursts where half of 
the primary are positrons; the large number of primary 
positrons enhances the positron annihilation flux enormously. 

It is interesting that, with the medium starburst fraction his- 
tory, the fiducial GeV emission from star-forming galaxies has 
a similar shape to the Fermi observed GeV background (Fig- 
ure [U. We find our fiducial background has a 1 - 30 GeV 
spectral slope of Fi 30 = 2.40, compare d to 2.41 for the ob- 
served background (lAbdo et al.ll2010al) . The predicted and 
observed spectrum only strongly diverge at energies below 
those observed with Fermi-LAiT, where the other backgrounds 
are measured with different instruments. However, this simi- 
larity of the GeV spectra is sensitive to the starburst fraction: 
for the low starburst fraction, the predicted 7-ray background 
from star-forming galaxies (Fi-so = 2.55) is softer than the 
observed spectrum; for high starburst fraction, the predicted 
7-ray background from star-forming galaxies (ri_3o = 2.27) is 
harder than the observed spectrum (see Figure [TO). Thus the 
apparent similarity may be coincidental. 

Summarizing our fiducial results, we find that starbursts and 
normal galaxies contribute roughly equally to the 1 GeV back- 
ground (— 17% each), with starbursts being more important at 
higher energies, but normal galaxies being much more impor- 
tant at MeV energies. 

However, there are a number of caveats with this predic- 
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tion. Most importantly, we have assumed that there are only 
two types of galaxies: those like the MiUcy Way, and those like 
M82. Essentially, we assume that the only kind of gas density 
evolution is the evolution of the starburst fraction. However, 
normal star-forming galaxies like the Milky Way had higher 
gas densities and star-formation rates (and therefore radiation 
fields) in the past. Very weak starbursts like the Galactic Cen- 
ter region may have hard TeV 7-ray spectra like M82 and 
NGC 253, but very low 7-ray luminosities because of winds. 
Very extreme starbursts like Arp 220 would probably be true 
proton calorimeters and be even more efficient at producing 
7-rays than M82 and NGC 253's core, although these are a 
small fraction of the cosmic star-formation rate. Submillime- 
ter galaxies have larger gas scale heights; if the CRs samples 
this gas evenly. Inverse Compton will be enhanced relative to 
bremsstrahlung. Finally, we have ignored IC emission from 
that have escaped normal galaxies, which will increase 
their MeV emission at low z. 

Comparison with previous predictions - As a check, we can 
compare our predictions to other recent predictions of the c os- 
mic GeV background. 

Using the high starb ur st f raction history, 
[Thompson. Ouataert. & WaxmanI (l2007l) predicted that 
~ 20% of the Fermi 1 GeV 7-ray background 10% 
of the EGRET GeV background) is from starbursts. We 
in fact predict ^3-9 times more GeV 7-rays with the 
T07 starburst history; with the fiducial starburst template, 
starbursts make up 58% of t he Fermi 1 GeV background. The 
differe nce arises because [Thompson. Ouataert. & WaxmanI 
((2007^ assumed much smaller values of 77 than we find, thus 
underpredic t ing th e 7-ray luminosity of M82 and NGC 253. 
iLacki et aLi ( 1201 lb noted this effect and argued that a high 
starburst fraction would lead to a > GeV starburst contri- 
bution of order ^ 50%, more in line with our predictions 
here. 

The pre-ferm/ calculation of iBhattacharya & Sreekumatj 
((2009) essentially scales the 7-ray luminosity of a galaxy to 
its TIR luminosity (roughly proportional to the star-formation 
rate), with different constants of proportionality for normal 
galaxies and starbursts, an approach comparable to ours. 
They find the total 7-ray background by integrating over in- 
frared luminosity functions. They find that only ^ 1% of 
the > 0.1 GeV background is from normal galaxies. This 
is probably b e cause they used the luminosity functions of 
iLa^ache et al.l (l2003h . in which normal galaxies make up 
only a small fraction of the cosmic star-formation rate at 
high z. (For comparison, the T07 starburst fraction gives 
us a '-^ 2% normal galaxy contribution to the 100 MeV 7- 
ray background.) They also used a hard E"^-^ spectrum for 
normal galaxies, although that only increases their predic- 
tions. They predict the starburst contribution to be < 6%; the 
small estimate likely arises because they derive (ultimately 
from equipartition argu ments about M82's radio luminosity 
from lAkyuz et al.lll991l with a volume scaling applied) that 
Lcev /iiR is only 30% higher for starbursts than for normal 
galaxies, whereas the Fermi observations indicate Lq^m jLx^ 
is abo ut 5-10 times higher 

Fields et alJ (I2OIO ) calculated the pionic 7-ray background 
from normal star-forming galaxies, relating the 7-ray emissiv- 
ity of a galaxy to its gas mass, which in turn is related to the 
star-formation rate of a galaxy through the Schmidt Law and 
empirical trends in the average radius of galaxies. Their vary- 
ing assumptions about what drives the evolution in the cosmic 
star-formation rate lead to predictions that normal galaxies 



are at least ^ 20% and possibly all of the unresolved 7-ray 
background. Because they used SFRmw = 1 M0 yr"', their 
normal galaxy contribution to the 7-ray background is twice 
as large as we would predict. After taking this in to account, 
the "density evolution" case of Fields et al. (2010) gives sim- 
ilar prediction of the normal galaxy contribution as us, while 
the "luminosity evolution" case with increased gas masses 
at high z predicts a higher normal galaxy contribution. The 
increased gas masses might make high-z normal galaxies as 
proton calorimetric as starbursts, in which case the luminos- 
ity evolution case predicts a similar contribution to the GeV 
background as the "starbursts" in our high starburst fraction 
m odel. 

iMakiva et all (1201 Ih predicted that - 5 % of the > 1 00 Me V 
7-ray background is from (mostly normal) star-forming 
galaxies, if the 7-ray luminosity per unit star-formation is the 
same as in the Milky Way (as we assume). This is some- 
what lower than our prediction of 6 - 17% of the 0.1-10 GeV 
background coming from normal galaxies. They also find 
that if gas evolution is taken into account, the contribution 
increases to ~ 10%. By using the Mitak a models of galaxy 
evolution (iNagashima & Y oshii 2004; N agashima et al.r2005l) 
and by approximating the starburst gala xy spectrum \yith the 
Milky Way spectrum plus a power law, i Makiva et al.l (1201 Ih 
also calculated the starburst contribution to GeV background. 
They found the starburst contribution was small, about 1 % of 
100 MeV background and a few percent of the 100 GeV back- 
ground. According to the Mitaka models, sta rbursts are only a 
few p ercent of the z w 1 star-formation rate (iNagashima et al.l 
120051) . which is of the same order as our low starburst frac- 
tion (HIO). By comparison, we find using the HIO starbust 
fraction that starbursts are ^ 0.6- 1% of the 100 MeV back- 
ground and 2-4% of the 1 GeV background, for the fiducial 
hi gh B models for M 8 2 and NGC 253's core. 

IStecker & VentersI (1201 1) performed a similar calcula- 
tion of the GeV ba ckground from star-forming galaxies to 
iMakiya et alJ (1201 ll) . although they used the HIO galaxy evo- 
lution models, and they assumed that starbursts had a soft. 
Milky Way-like spectrum. They found that normal star- 
forming galaxies were anywhere between ^ 10% and all of 
the GeV background depending on how the 7-ray luminosity 
scales with star-formation rate and/or gas mass. Since they 
include the effects of gas evolution for normal galaxies, it is 
natural that their results would be higher than ours. They find 
a starburst contribution peaking at ^ 1% to the GeV back- 
ground, which is somewhat lower than our results for an HIO 
star-formation rate, and dropping off more quickly at high en- 
ergies than we predict because of the soft 7-ray spectrum they 
used. 

Overall, given the differing methods of other calculations 
and some large uncertainties in the input parameters (particu- 
larly the starburst fraction), our results are broadly consistent 
with previous predicti ons. 

We also note that IChakraborty & Fieldsl (1201 2l) . concur- 
rently to our paper, have recently presented a calculation of 
the 7-ray background from IC and pionic emission from nor- 
mal galaxies. 

5.4. Constraints from spectral shape 

The star-forming galaxy contribution to the GeV 7-ray 
background remains unclear, despite the detection of several 
star-forming galaxies (both normal and starburst) by Fermi. 
In addition, our predictions rely on assumptions about the 
star-formation rate of the model galaxies, and the Milky Way 
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Fig. 1 1 . — Predicted 7-ray background from normal galaxies (solid) using 
our medium starburst fraction with the h = 1000 pc (grey) and the h = 2000 pc 
(black) models, as compared with the observed 7-ray background (blue). 
Also shown are the spectrum normalized to fit the observed background 
(dashed). Another large component at energies below 100 MeV and a compo- 
nent at energies above 10 GeV are needed to get the shape of the background 
correct. 

needs an additional scaling factor to match its known 7-ray 
luminosity. However, star-forming galaxies obviously con- 
tribute < 100% of the GeV background. Assuming we know 
the SED of star-forming galaxies from MeV to GeV energy 
ranges, we can therefore set an upper limit on the star-forming 
galaxy contribution to the MeV background by normalizing 
the SED so it touches but never exceeds the GeV background. 

Table |4] lists the maximum fractions of the observed MeV 
background using different models normalized to the ob- 
served GeV background at finonn- 

The shape of the normal galaxy 7-ray SED is incompatible 
with normal galaxies being the majority of the MeV back- 
ground, even with the increased leptonic emission from the 
CMB at high z- The shape of the spectrum constrains nor- 
mal galaxies to be a half or less of the 100 MeV background, 
< 25% of the 30 MeV background, < 10-20% of the 10 MeV 
background, and < 10% of the 1 MeV background. Note that 
even these values require the normal galaxy contribution to 
be scaled up by a factor ~ 6-25. The discrepancy in shape is 
shown in Figure fTTl there is a large hump of additional energy 
near 10 MeV that normal galaxies cannot explain. 

For starburst galaxies, the MeV contribution is typically 
even smaller, although the exact fraction of the MeV back- 
ground depends on the magnetic field strength. The contribu- 
tion from starbursts does not need to be multiplied by much 
(factor < 2) to account for the GeV background, as long as the 
starburst fraction is high; with the medium starburst fraction, 
the emission must be scaled up by a factor ^2-1 to reach 
the GeV background. In high B (hadronic) models, starbursts 
are less than 3-5% percent of the 30 MeV background and 
only < 2% of the 10 MeV background (black lines in Fig- 
ure [T2]l. Only a fraction of one percent of the 1 MeV back- 
ground could be nonthermal emission from high B starburst 
galaxies. Low B (leptonic) models have larger allowed MeV 
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Fig. 12. — Predicted 7-ray background from starbursts (solid) using our 
medium starburst fraction with the B = 300 fiG (black) and B = 50 (grey) 
best-fit models of M82. The thick fines assume that primary are all elec- 
trons, while the thin lines are if the primary are half positrons, half elec- 
trons. The observed 7-ray background is shown in blue. The dashed lines are 
tfie predicted spectra normalized to fit the observed background. Below 10 
GeV, another very large component is needed to get the shape of the back- 
ground correct. 

fractions, but still make up a small minority of the MeV back- 
ground: < 15-25% of the 30 MeV background, < 10-20% 
of the 10 MeV background, and < 3-8% of the 1 MeV back- 
ground (grey lines in Figure [12). Including inflight positron 
annihilation radiation from primary positrons increases the 1 
MeV contribution, but not enough to be a major contributor 
of the background. As Figure [12] shows, starbursts cannot be 
the main source of the MeV background either 

We therefore conclude that star-forming galaxies cannot 
make up a significant portion of the MeV background, even 
if the overall 7-ray production efficiency is much higher than 
we suppose. Either (1) some other source of MeV emitters 
is needed, (2) some unknown emission process make star- 
forming galaxies much brighter in MeV emission than we 
expect (and much brighter per unit star-formation than the 
Milky Way appears to be), or (3) the claimed cosmic MeV 
background is incorrect. 

6. CONCLUSIONS 

The origin of the MeV background remains mysterious, as- 
suming the claimed intensity of the background is correct. 
Star-forming galaxies are an attractive source for the GeV 
background because they are ubiquitous, the physics is at least 
partly understood (especially in the Milky Way), and several 
have now been detected by Fermi and TeV telescopes. How- 
ever, most of their GeV emission is likely pionic, and there- 
fore rapidly falls off below 100 MeV; hence this process 
cannot be responsible for the MeV background. We have con- 
sidered possible processes that can provide MeV emission, 
including (1) leptonic Inverse Compton and bremsstrahlung 
emission produced by CR in star-forming galaxies, (2) 
7-ray line and positron annihilation radiation, and (3) inter- 
galactic cascades that shift energy from the TeV band to the 
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MeV-GeV band. The IC is expected to be most important in 
normal, quiescent galaxies like the Milky Way, because CR 
protons escape more easily and thus primary CR electrons 
are more important; furthermore, the growing strength of the 
CMB radiation field will strengthen IC losses in high z normal 
galaxies. Intergalactic cascades, conversely, will be more im- 
portant for starburst galaxies which are observed to have hard 
spectra extending to the TeV. 

We have modeled the Milky Way and the starbursts M82 
and NGC 253 placed at various redshifts. From these mod- 
els, we have assembled the expected spectrum of the 7-ray 
background, taking into account redshift and intergalactic cas- 
cades. Our primary results are: 

• The MeV emission of normal galaxies is greatly en- 
hanced by the increased energy density of the CMB 
at high z- This increases Inverse Compton losses at 
the expense of other radiation and cooling mechanisms 
like synchrotron. As a result, the pion bump of normal 
galaxies is "filled in" at high z (Figure[3]l. 

• Starburst galaxies have relatively low MeV emission 
compared to their strong GeV emission. The relatively 
strong proton losses imply a strong pion bump, since 
there is more power in CR protons than in CR e^. Cas- 
cade emission increases the MeV emission by a factor 
< 2, nowhere near enough to fill in the pion bump. This 
conclusion can be altered in "low B" models of star- 
bursts, in which the GeV emission comes mostly from 
Inverse Compton emission (Figures|5]and|8}. 

• We find with our "medium" starburst fraction (15% of 
star-formation at all z) that star-forming galaxies are 
^ 1/3 of the GeV 7-ray background, with roughly 
equal contributions from normal galaxies and starbursts 
(Figure [U. However the effective starburst fraction is 
highly uncertain and depends on gas evolution in nor- 
mal galaxies. As the fraction of cosmic star formation 
in starbursts approaches unity, the observed GeV 7-ray 
background is produced or even exceeded (Figure [TO). 
Starbursts dominate the 7-ray background at high ener- 
gies, while normal galaxies dominate at low energies. 

• In our fiducial model, star-forming galaxies make up 
only 1-2% of the 1 MeV background and 4% 
of the 10 MeV background (Figure [T]!. Much of this 
emission is IC from normal galaxies, but there are also 
important contributions from 7-ray lines from CCSNe, 
positron annihilation, and ^^Al decay. 

• To account for the uncertainties in 7-ray luminos- 
ity normalization, we derive more robust limits on 
the MeV contribution from the spectral shape of star- 
forming galaxies. According to these constraints, nor- 
mal galaxies could make at most half of the 100 MeV 
background, 25% of the 30 MeV background, and 20% 
of the 10 MeV background (Figure fTTT i. However, 
reaching these limits requires the normal galaxy back- 
ground to be multiplied by a factor ^10. Similarly, 
high B starbursts could make up at most ~ 15% of the 
100 MeV background, 5% of the 30 MeV background, 
and 2% of the 10 MeV background. If the starburst 
emission is leptonic (B w 50 /iG), these fractions could 
be higher, with ^ 10-20% of the 10 MeV background 
allowed to be from starbursts (Figure [12) • 



In short, while star-forming galaxies can (and according 
to our models, do) provide a substantial fraction of the GeV 
7-ray background, they are a minor contributor of the MeV 
background. The reason why the MeV background cannot 
be explained by star-forming galaxies is simply that the 7- 
ray background is very steep: a lot more power is required 
at MeV energies than at GeV energies. Thus, even if star- 
forming galaxies do have enough luminosity at GeV energies, 
and even if the spectrum is flat at MeV energies (as is the 
case for normal galaxies), this is insufficient for the observed 
background. 

Thus, there must be some other explanation for the MeV 
background. The simplest explanation would be some other 
source that peaks at MeV energies. But this just again raises 
the mystery of why the MeV-GeV 7-ray background appears 
to be a featureless power law if there are different sourc es re- 
sponsible at different energies (IStecker & Ventersll201 lb . It is 
possible this is sheer coincidence. Indeed, something like this 
happens at GeV energies in our fiducial model: normal galax- 
ies with soft spectra are most important at energies below a 
few GeV, while starbursts with hard spectra are most impor- 
tant above a few GeV; the sum is a nearly featureless power 
law with r = 2.4, as observed. Likewise, current blazar mod- 
els of the 7-ray background invoke a coincidence with soft 
Flat Spectrum Radio Quasars at low energies and hard BL 
Lacs at high energies (e.g.. Venters & Pavlidou 2011). We 
should also keep in mind that the featureless power law in 
Figure [T]is actually for the currently unresolved background. 
The distinction between unresolved and resolved sources is 
an experimental one (see the huge sensitivity gap in Figure [13] 
between COMPTEL and ferm/-LAT), not necessarily one 
with fundamental astrophysical significance. The resolution 
of sources is much better at GeV energies than MeV energies, 
and the resolved 7-ray blazars would add a large component 
on top of the unresolved background. Finally, the reported 
MeV background could simply be incorrect (due to detector 
backgrounds) or Galactic in origin. 

Our conclusions can be evaded if there are extremely high 
amounts of low energy CRs below the pion production thresh- 
old, or extremely large contributions from discrete MeV 
sources. IC emission from escaping in the halos around 
normal galaxies can increase the MeV emission, but since the 
IC emission is a broad continuum extending to GeV ener- 
gies, the spectral shape constraints imply that this emission 
does not make up all of the MeV background. Low energy 
spikes of CR electrons can enhance the bremsstrahlung and 
IC emission, but can be constrained with synchrotron radio. 
An extreme enhancement of MeV nuclei might even increase 
the nuclear line emission enough to alter the SED, although 
the ionization rate from these nuclei should also be extreme. 
Discrete sources are not expected to be dominant, based on 
observations of the Milky Way, so a new significant MeV 
source population would have to scale non-linearly with star- 
formation rate. 

There are several indirect tests of our predictions. The 
relatively strong MeV contribution of normal galaxies arises 
from increased Inverse Compton losses on the CMB, which 
comes at the expense of other losses - particularly ra- 
dio. Thus at high z, it has long been expected that 
galaxies should be faint compar ed to expectati o ns from 
their infrared luminosities (Carill i & YunI Il999t iMurphvl 
120091) . In contrast, observations indicate that galaxies lie 
on the empir ical z = infrared-radio correlation to relatively 
high z (e.g., lAppleton etaT][2004t iGam & Alexaiide^ [20091 
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Fig. 13.— Fiducial spectra of M82 and NGC 253's core (black solid) and 
M31 (black dashed) compared to existing (grey) and proposed (light grey) 
7-ray instruments. Even with a 10^ increase in MeV sensitivity over COMP- 
TEL, detection of these galaxies will be difficult. The sensitivity curves 
assume exposure times of 10* seconds and that A£ = E, with data taken 
from IKa nbach et al. (2004), Tajima et al. 1 2005). iBloser et alj <200^ . and 
IGreiner et al.i t2011i) . The proposed MEGA instrument would have a simi- 
lar sensitivity as GRIPS. 



ISargent et al]|20 T0^,'Mao et al."201 A or even are radio-bright 
re.gjKovacser al. 2006; Vlahakis et al. 2007; Murphy 2009; 
llvisonetalj|201 0; Micha towski et al.ll201Q) . At high z, these 
observations probe galaxies with high star-formation rates and 
possibly more starburst like, where the effects of the CMB are 
expected to be weak. However, new instruments like JVLA 
and eventually SKA may probe Milky Way-like galaxies, and 
determine the effects of the stronger CMB. 

Further measurements extending the hard X-rays and the 
high energy 7-rays towards the MeV range may prove con- 
straining for the starbursts. ferm/-LAT should ultimately 
reach down to 30 MeV. At these energies the pion bump 
should fall off in the high B models, but the leptonic emission 
peaks in the low B models. Thus /rrm;-LAT can help con- 
strain B in the nearest starburst s. Further Fermi-LAT d ata will 
resolve out additional blazars (lAbazaiian et al.ll201 Ih . which 
would constrain the true contribution of star-forming galax- 
ies to the GeV background and limit their contribution to the 
MeV background. In the hard X-ray regime, measurements 
of the 10 - 100 keV luminosity of starbursts would limit their 
IC emission, helping constrain B, and synchrotron emission, 
which can be significant to ^ 1 MeV in some of our models. 

Better MeV instruments may also help us understand the 
emission mechanisms of 7-rays in this energy range within 
the Milky Way, particularly the contribution of discrete 
sources. As noted in section lTSl both GalProp models and our 
own models (especially the h = 1000 pc model) fail to predict 
all of the emission around ^ 10 MeV (Figure|2]l. The discrep- 
ancy is poorly understood, and implies there is something we 
do not know about galactic MeV emission. The discrepancy 
by itself cannot explain the MeV background - the Milky Way 
is less luminous at 10 MeV than at 1 GeV, unlike the back- 



ground - but it may increase the predicted background by a 
factor of a few. 

However, direct measurements of the 1-30 MeV emission 
of star-forming galaxies beyond the Milky Way and possi- 
bly the Magellanic Clouds seem far off. In Figure [13] we 
plot the fiducial models of M82 and of NGC 253's core com- 
pared to the sensitivities of current (grey) and possible future 
(light grey) 7-ray instruments. In X-rays and GeV-TeV 7- 
rays, the highly sensitive instruments we have now allow de- 
tection of these starbursts. However, COMPTEL, the best in- 
strument over much of the MeV band, is about three orders 
of magnitude behind ferm;-LAT. Nor wo uld proposed next- 
gener ation instru ments like MEGA (e .g., IB loser et al.l 120021 
I2OO6I) or GRIPS (IGreiner et al.l |20 111) be able to detect the 
nonthermal emission from M82 or NGC 253. But even a 
MeV instrument a s sensitive as FermZ-LAT, such as ACT 
jMilne et alj|2002l; iBoggs 2006), might not detect M82 and 
NGC 253 from 1-10 MeV: the pion bump in the fiducial 
model means that these galaxies are much fainter in the MeV 
band than in the GeV band. We also consider M3 1 by scal- 
ing the Milky Way fiducial model to the star-formation rate 
of M31 ( 1 Mf7) yr~'; I Williams 2003) and using a distance of 
780 kpc (IStanek & Gai-navich,1998) . While M31 has a softer 
MeV-GeV spectrum, it also has a much lower 7-ray spectrum 
to begin with, so it too is difficult to detect at MeV energies. 
Thus, vastly more powerful MeV telescopes are required to 
directly observe the MeV emission of individual star-forming 
galaxies; such improvem ents may be feasible with new 7-ray 
optics technologies, as in lHabs et alJ (l2012l) . 

The origin of the MeV (and GeV) background can be tested 
by th e anisotropics in the background (Zhang & Beacon^ 
I2OO4I) . Normal gala xies are expected to ha ve a very weak 
anisotropy signal ( Ando & Pavlidoul 120091), wh ile blazars 
have strong anisotropy signals (lAndo et al.l 120071) . The sig- 
nal for starbursts has not yet been calculated, but is probably 
stronger than normal galaxies; anisotropics in the submillime- 
ter associated vyith su bmillimeter galaxies have been detected 
dAmblard et al.ll20Tlb . At GeV energies, the strength of the 
star-forming background is highly sensitive to the starburst 
fraction; more proton-calorimetric galaxies are brighter in 7- 
rays. Thus, there might be strong anisotropics in the GeV con- 
tribution from star-forming galaxies if starbursts dominate. In 
the MeV, however, normal galaxies are expected to dominate, 
and these have very weak anisotropy signals. A future MeV 
experiment could therefore search for the star-forming galaxy 
contribution by looking for a uniformly distributed compo- 
nent to the MeV background. Further studies of the starburst 
contribution to background anisotropics, and a better treat- 
ment of galaxy evolution, is necessary to evaluate this effect. 

In short, the sheer power of the MeV background as re- 
ported requires many luminous sources, and star-forming 
galaxies simply are not bright enough to be its origin. How 
the MeV and GeV backgrounds join up so seamlessly is un- 
clear, since star-forming galaxies seem to make up at least 
^ 20% and plausibly a majority of the GeV background. Fu- 
ture MeV instruments must be improved, possibly by many 
orders of magnitude, to solve the riddle of its source. 

Note - As we w ere completing this paper, w e became 
aware of the paper bv lChakrabortv & Fieldsl(l2012l) that calcu- 
lates the contribution of Inverse Compton 7-rays from normal 
galaxies to the 7-ray background. 
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APPENDIX 



General approach - Our approach to one- zone models is ultimately ba sed on the extensive discussion in lTorresI (l2004l) . and is 
described in detail in lLacki et al.l (1201 Ol) and lLacki & Thompson! (1201 Obi) . 
The population of CRs in galaxies is governed by the diffusion-loss equation: 

-DV^N(E) ^^-^ WEmE)] - QiE) = (1) 
t(E) dE at 

where D is the spatial diffusion constant, b(E) is the sum of all of the energy loss rates, Q{E) is the CR inje ction rate for CRs a t 
the energy E, t(E) includes all catastrophic and escape losses, and A^(£') is the CR energy spectrum (see also lStrong et al.ll2007b . 
Our steady-state one-zone models assume the CR population has no time or spatial dependence within the modeled region, so 
that the diffusion-loss equation reduces to the leaky box equation: 

N{E) d 

-^-—[b(E)N(E)]-Q(E) = 0. (2) 

This may be solved by the Green's function method described in lTorresI (|2004|) . 

CR nuclei heavier than hydrogen - In order to estimate the strengths of the CR de-excitation nuclear lines, we treat primary 
helium-4, carbon-12, and oxygen-16 ions in addition to protons, electrons, and positrons. However, secondary nuclei are not 
included. 

We normaUze these nuclei by the ratio of their rigidity spectra with CR protons' rigidity spectrum, comparing at equal rigidities. 
For our momentum power law spectra dQ/dq = Cq~'', it can be shown that 

Ci _ {dQ/dR)i (Z;y-' 



Cp (dQ/dR)H \Zp 

where R = qc/{Ze) is the rigidity, and Z is the electric charge o f the nucl eus (Z„ = 1). The ratio of the helium -4 and hydrogen 
rigidity spectra is assumed to be 8, as observed in the Milky Way dWebber & Lezni akll974tlWebber et al.ll987l) . The abundances 
of carbon and oxygen are then normalized to helium-4 using the values in Me ver et al. 

Heavier nuclei can interact with gas in the ISM and not only produce pio ns, but undergo nuc lear reactions as well. To calculate 
these lifetimes, we use a modified version of the inelastic cross sections of iLetaw et al.l ( Il983h . The cross section for interactions 
of CR nuclei and ISM protons is 

af^^^{E) = A5A^(P^[\+Qm6 sin(5.3-2.631nAcR)][l -0.62e-^/2"'"^'=^sin(10.9(/i:/MeVr°-^'*)] mb (4) 

for CRs with an atomic mass Aqr and kinetic energies K above 10 MeV. For interactions between CR nuclei and ISM helium 
nuclei, we replace the term with [4'^'^^ +^cr'' ~ 1]^ to get a cross section a-^^{E), by analogy with the Glauber rule. The 
lifetime of a CR nucleus to nuclear reactions is then 

fnuclear = [«el(£)«H + <d(£)nHe)/3c]-' , (5) 

where /3c is the speed of the CR. 

Pionic emission from nuclei - While pionic cross sections for proton-proton collisions are given in iKamae et al.1 (l2006l) . no 
cross sections are given for proton-helium or helium-helium collisions in that reference. Since the helium is expected t o be a 
minor contributor to the pionic emission 10%), we scale the proton-proton cross sections using the lOrth & Buffingtonl(ll976h 
approximation: 

af (£z) = (£z/Acr)[Ac/r +Ai/m- 1]' (6) 

which is applied to the differential cross sections for production of pionic 7-rays, e^, and neutrinos.'^ Roughly speaking, this 
rule takes into account that nuclei are not transparent to hadronic interactions, with some of the nucleons being shadowed by the 
others. 

Nuclear lines - Nuclei colliding with each other can excite nuclei; the de-excitations then generate 7-ray lines in the rest-frame 
of the emitting nucleus. There are three different cases: a CR nucleus can hit an ISM proton and become excited, a CR proton can 
hit an ISM nucleus and excite it, and a CR nucleus can penetrate a dust grain. Excited CR nuclei will generate a broad continuum 
as their Doppler shift smears out the 7-ray line, whereas excited ISM nuclei will generate a narrow line. We neglect the dust grain 
case, but it would generate an even narrow er line as energy losses in the dust grain dampened the nuclear recoil. 

FoUowing lMeneguzzi & ReevesI (Il975h . we calculate the line emissivity of CR nuclei on ISM p/a collisions as 



J dE' dE"^(E')ca,i{E')-^5(E"-EL), (7) 



dE 2^"*7 " 7w, ~ dE'- '7' 

k,i '^min 

■ The lOrth & Buffingtor] il97(i) approximation is also used by the GALPROP code for simulations of the Milky Way (Moskalenko & Strond 19981) . 
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where the emissivity is in units of ph cm"'' s~' ergs"' (contrast with lMeneguzzi & ReevesI (Il975h . where it is photons per ISM 
hydrogen atom per steradian), E'^^^ = £'/(7(l + and £,"3^ = £'/(7(l -/3)). 
For the hne emissivity of CR p/a on ISM nuclei, we use 



dE ^ 2AEl 

k.i 



J dE'ai,(E')cfi\^{E') X {\E-El\ < A£J. (8) 



For simpHcity, we assume the nuclear line has a line half-width of AE^, which we take to be 40 keV. Since we are interested in 
the cosmic background, redshift will sm ooth out narrow feature s in any case, so the exact line width is not that important. 

We use the line cross sections from (iKozlovsky et alJ [20021) . Included lines are 429 and 478 keV lines from '^He (broad 
components only since the projectile and target are both the same), 4.438 MeV and 15.1 MeV lines from '^C (where the 4.438 
MeV line can be created in collisions with '^N and '^O as well because '^C is created in the collision) , and 6.129 MeV, 6 .916 
MeV, and 7. 115 MeV Hnes from "'O. The 4.438 MeV '^C line is generally expected to b e the strongest (iRamatv et al.|[T979h . 

Bremsstrahlung - In this paper, our models use the bremsstrahlung cross sections from lStrong et al.l(l2000ab . An electron with 
kinetic energy interacts with an atom of number Z and atomic mass A^, produces a photon with energy krUgC^, and leaves with 
a kinetic energy Ke - k nieC^ > 0. We inc lude the effects of hydrogen and helium in the ISM; metals have a negligible effect on 
the total cross section. IStrong et al.l (l200 0a) consider three regimes for the electron kinetic energy: 10 keV - 70 keV, 70 keV - 2 
MeV, and more than 2 MeV. They also consider the Fano-Sauter limit when E^ « K^, because the other approximations to the 
bremsstrahlung cross section drop off too quickly. We use the Fano-Sauter limit when it gives a bigger cross section than the 
other approximations, although we find in practice this only happens when km^c^ is very nearly K^. 

The bremsstrahlung loss time is then calculated by directly integrating up the cross sections: 

(ISchlickeiseiil2002h . 

Positron annihilation - Positrons have a small but significant probability of annihilating in-flight, before thermalizing. The 
positron annihilation cross section is given as 



((r^ + 4r-8)log ^+^ -(r+4)v/7(7^) (10) 
t^(t-4) \ 4 / 



where r = 2(1 + 7) fe-g - lBeacom & Yukselll2006tlVietrill2008h . Then the in-flight annihilation time is (neCTann/^ec)"', where ne is 
the number density of both free and bound electrons. 
The inflight annihilation spectrum is then 



7+ + 1 - £7 £7 / V £7 7+ + 1 ~ £7 / V ^7 7+ + 1 



(11) 



where e^, = E^/inieC^), p+ = y^j^-l, and 7+ > (l-2e^ + 2e^^)/(2e^- 1) (e.g.,'Ahar onianll2004l) . 

Positrons are more likely to cool by ionization and annihilate after thermalizing than at relativistic speeds. Positrons are 
conveyed towards rest in energy space at a "speed" \dE/dt\. The rate at which positrons annihilate at low energies in equilibrium 
is therefore Q^" = dN /dE x \dE /dt\ x fsm-v(E), where /sui-v(£^) is the probability that positrons survive cooling to rest without 
being catastrophically lost (through inflight annihilation or escape). We take our lowest energy bin (with kinetic energy Kq = 
100 keV) and assume /sui-v(^o) = 1 to calculate the positron annihilation rate as 

dE 

Note that this is not the rate at which positrons are injected at all energies: positrons can escape the system through diffusive and 
advective losses at high energy and avoid annihilation, or annihilate in-flight. 

Annihilation at rest can either occur through direct annihilation or through positronium formation. Direct annihilation leads to 
two 511 keV line photons. Positronium formation, which occurs /ps sa 90-95% of the time in the Milky Way, also leads to two 
511 keV line photons about 1/4 of the time. We arbitrarily assume a line half-width of AE = 5 keV (which will be smeared out 
for galaxies at different z anyway) to calculate the 511 keV annihilation line from CR positrons as: 

dQ ^ (l-3M/4)Qr 

dE 2AE ^ ' 

We assume /ps = 0.9 for this work. Positronium annihilates into three photons with a continuum distribution of ener gies b elow 
511 keV 3/4 of the time. The positronium continuum spectrum is given in Ore & Powell ( 1949) and Guesso um et al.. (iI991h : 



dQ _ 6 

dE (tt^ - 9)m^c^ 



£(!-£) , 2(l-e) , oli-fT, n 

7S — ^ + 2— — -3log(l-e) + — 

{2-ey (2-e)-* e 



3Gr"/ps 



(14) 



where e = E /(nieC^) and we have already taken into account that each annihilation produces three photons. 



MeV BACKGROUND FROM STAR-FORMING GALAXIES 



25 



77 Absorption - At TeV ener gies, 77 absor ption process ma y become important. This p rocess converts 7-rays into e^, which 
can radiate in Inverse Compton (Ilnoudl201 Ibl) or synchrotron (iLacki & Thompsonll2010bl ). The interior 7-ray photon density is 



calculated using the uniform slab model, where the absorbing and emitting regions are cospatial: 

Nj = -^^[1 -exp(-T^^(/i))] (15) 

where Qj is the rate at which 7-rays are injected per unit volume, h is the midplane-to-edge scale height, and T-^^{h) is the 
midplane-to-edge optical depth to 77 absorption. To calculat e T^^jh) = h J n(e) a~/~,{e,E^)de for a radiation field n(e), we use the 
approxi mation for the 77 cross sections a-y-y(e,Ey) given in lAharonianl (f2004i) . The spectrum of these pairs is then calculated 



using the lAharonian et all ( Il983h source functions. 

For the Earth-observed spectrum, we again use a uniform slab model, but substitute the galaxy radius R for h, since the Milky 
Way, M82, and NGC 253 are all viewed nearly edge-on: 

F® = -exp(-r^^(7?))]/T^^(7?). (16) 

where i^™^*'^ is the unabsorbed 7-ray flux. 

Finally when constructing the background, we assume that a galaxy is typically observed face-on. This means that the sightline 
within the disk will typically have length of order ^ h. However, the radiation field of a galaxy extends out to a distance ^ R> h. 
Thus 7-rays in the halo of the galaxy/starburst may be absorbed by the 77 process. The pair in the halo will not radiate in 
bremsstrahlung; though they will radiate IC and possibly synchrotron, depending on how far out the magnetic field extends. To 
be conservative we do not include any emission from these halo e*. Instead, we apply the uniform slab model with sightline h to 
the galaxy itself, and then treat the halo as a foreground screen of length R: 

F^'^^He,) = F---- 1 -^^P(-^y(^)) exp(-r,,(j;)), (17) 
where we assume the radiation field is constant out to a radius R. 
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2.7 


3.7 


8.3 


15 


22 


28 


15 






TO? 


0.18 


0.75 


3.4 


13 


15 


20 


44 


76 


110 


140 


57 






HIO 


0.035 


0.14 


0.16 


0.58 


0.59 


0.79 


1.7 


2.9 


4.2 


5.4 


2.5 






15% 


0.14 


0.84 


0.81 


2.8 


2.9 


3.9 


8.4 


15 


22 


28 


15 






TO? 


0.85 


3.7 


4.3 


15 


16 


21 


45 


77 


110 


140 


59 


B = 150 IJ.G 


e~ 


HIO 


0.004 


0.017 


0.071 


0.26 


0.3 


0.43 


1 


1.9 


2.9 


4.4 


2.5 






15% 


0.019 


0.082 


0.35 


1.3 


1.4 


2.1 


5.1 


9.6 


15 


23 


15 






T07 


0.11 


0.44 


1.9 


7 


7.9 


11 


27 


50 


77 


120 


59 






HIO 


0.018 


0.075 


0.089 


0.3 


0.32 


0.44 


1 


1.9 


2.9 


4.4 


2.5 






15% 


0.072 


0.44 


0.44 


1.5 


1.5 


2.2 


5.2 


9.6 


15 


23 


15 
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Model 


Prim. 


/sb'' 










Predicted fraction of 7-ray background 
















100 keV 


300 keV 


1 MeV 


3 MeV 


10 MeV 


30 MeV 


100 MeV 


300 MeV 


1 GeV 


10 GeV 


100 GeV 










% 


% 


% 


% 


% 


% 


% 


% 


% 


% 


% 








T07 


0.45 


1.9 


2.3 


7.9 


8.4 


12 


28 


50 


77 


120 


59 


B 


= 200 AiG 




HIO 


0.0035 


0.016 


0.073 


0.31 


0.39 


0.63 


1.7 


3 


4.3 


4.7 


1.9 








15% 


0.016 


0.076 


0.36 


1.5 


1.9 


3.1 


8.2 


15 


22 


25 


12 








T07 


0.09 


0.42 


1.9 


8.1 


10 


17 


44 


80 


110 


120 


44 






e± 


HIO 


0.037 


0.13 


0.12 


0.38 




U.DO 


1 7 


J 


4.3 


4.7 


1.9 








15% 


0.16 


0.75 


0.63 


1.9 


2 1 


3 2 


S 'K 

O.J 




22 


25 


12 








T07 


0.93 


3.3 


3.3 


10 


1 \ 


17 


45 


81 


110 


120 


45 


B 


= 250 /iG 


e' 


HIO 


0.01 


0.024 


0.068 


0.22 


27 


48 


1 4 


2 8 


4 


4.1 


1.5 








15% 


0.051 


0.12 


0.35 


1.1 


1.3 


2.3 


7.1 


14 


21 


22 


9.8 








T07 


0.26 


0.64 


1.8 


5.7 


7.3 


13 


39 


73 


100 


110 


36 






e± 


HIO 


0.009 


0.047 


0.042 


0.12 


u. 1 J 


90 
U.Zz 


u.o / 


1 'X 


2 


2.1 


1 








15% 


0.044 


0.29 


0.23 


0.64 


71 


1 2 


J.O 


/ .J 


11 


12 


6.8 








T07 


0.23 


1.2 


1.1 


3.1 


3 5 


5 9 


18 


35 


51 


53 


22 


B 


= 300 tiG 


e" 


HIO 


0.0022 


0.0099 


0.04 


0.16 


21 


39 


1 2 


2 4 


3.8 


5.4 


2.6 








15% 


0.01 


0.049 


0.2 


0.76 


1 


1.8 


5.7 


12 


19 


28 


16 








T07 


0.057 


0.26 


1.1 


4.2 


5.7 


10 


31 


63 


100 


140 


61 






„± 


HIO 


0.01 1 


0.047 


0.052 




0.23 


0.4 


1.2 


2.4 


3.8 


C /I 

5.4 


2.6 








15% 


0.044 


0.27 


0.26 


0.88 


1.1 


1.9 


5.8 


12 


19 


28 


16 








T07 


0.27 


1.2 


1.4 


4.8 


6.1 


11 


32 


64 


100 


140 


61 


B 


= 400 /iG 


e' 


HIO 


0.0016 


0.0074 


0.027 


0.1 


0.15 


0.33 


1.2 


2.6 


4.4 


7 


3.6 








15% 


0.0075 


0.037 


0.14 


0.47 


0.72 


1.5 


5.7 


13 


22 


37 


23 








T07 


0.042 


0.2 


0.71 


2.7 


4.1 


8.7 


32 


69 


120 


180 


85 








HIO 


0.004 


0.017 


0.03 


0.11 


0.16 


0.33 


1.2 


2.6 


4.4 


7 


3.6 








15% 


0.017 


0.098 


0.15 


0.51 


0.74 


1.5 


5.7 


13 


22 


37 


23 








T07 


0.1 


0.45 


0.79 


2.8 


4.2 


8.8 


32 


69 


120 


180 


85 



Denotes whether the primary electrons and positrons in the model are all electrons or half electrons and half positrons (e ). 

Starburst fraction used in the model. HIO: our low starburst fraction case, with the starburst fraction in the lHopkins et alj|201Ct) semianalytic models; 15%: our fiducial medium 
starburst fraction case, where 15% of the cosmic star formation rate at all z is in starbursts; T07: our high starburst fraction case, using the values in lThompson. Ouataert. & WaxmanI 
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TABLE 4 

Constraints From Spectral Shape 



Model 


Prim. 


/sb" 


TP C 
'^norm 


Scale"* 








Maximum fraction of 7-ray background^ 




















100 keV 


300 keV 


IMeV 


3MeV 


lOMeV 


30MeV 


100 MeV 300 MeV 


1 GeV 


10 GeV 


100 GeV 








GeV 




% 


% 


% 


% 


% 


% 


% 


% 


% 


% 


% 




Normal Galaxies: Milky Way 


h = 1000 pc 


e 


HIO 


1.1 


7 


0.9 


2.8 


8.5 


17 


15 


19 


49 


83 


100 


68 


20 




15% 


1.1 


8 


0.9 


2.8 


8.5 


17 


15 


19 


49 


83 


100 


68 


20 






T07 


1.5 


25 


0.41 


1.6 


6.6 


11 


10 


13 


36 


74 


98 


77 


35 


h = 2000 pc 


e~ 


HIO 


1.1 


6 


1.4 


3.9 


12 


27 


21 


24 


52 


83 


100 


69 


20 




15% 


1.1 


6.8 


1.4 


3.9 


12 


27 


21 


24 


52 


83 


100 


69 


20 






T07 


1.7 


21 


0.66 


2.2 


8.4 


18 


15 


18 


40 


74 


97 


82 


36 
















,Stai 


-bursts: MX2 
















— n 7t 

-D = 50 flG 




HIO 


A 

4 


o^; 
ZD 


A AA 

U.44 


1.0 


0.0 


00 
ZZ 


OA 

2U 


O/t 

24 


A1 

4i 


65 


00 
00 


A1 

yi 


io 




Vj/o 




4.5 


A A 

U.4 


l.D 


0.1 


1 

Zl 


1 

ly 


00 

22 


A 1 
41 


OJ 


50 


00 
y3 


1 

21 






lU/ 




u.yo 


A AA 
U.44 


1.0 


0.0 


00 


01 
21 


O/l 
24 


A'X 
4J 


00 


QQ 

00 


01 
yi 


1 A 
10 






HIO 


4 


26 


1.4 


5.5 


7.7 


25 


21 


24 


43 


65 


88 


91 


18 






15% 


4.3 


4.8 


1.1 


5.9 


7.2 


23 


20 


23 


41 


63 


86 


93 


21 






T07 


3.8 


0.98 


1 .3 


5.3 


7.8 


oc 
25 


00 
22 


O/l 

24 


A A 

44 


66 


88 


01 
91 


1 A 

16 


-D = lUU fJXJ 




HIO 


3.5 


39 


0.26 


1 


A O 

4.Z 


15 


1 

16 


1 
21 


A A 

44 


OA 
/O 


AO 

93 


Ac 

95 


00 

33 






1 

3.0 


7.3 


0.23 


A AO 

0.93 


A 

3.9 


1 A 

14 


15 


1 n 
19 


A 1 

41 


68 


92 


OA 

96 


Af\ 

4u 






lU / 




1 c 
l.J 


A O/^ 
U.ZD 


1 
1 


A O 
4.Z 


1 c 
ID 


1 A 
10 


1 
21 


A A 
44 


"71 
/ 1 


00 
yj 


oc 

yj 


OA 

3U 




± 

e 


HIO 


3.5 


39 


0.92 


3.8 


5 


17 


17 


22 


44 


71 


93 


95 


33 






15% 


3.8 


13 


0.7 


4.1 


4.7 


16 


15 


20 


42 


68 


92 


96 


40 






T07 


3.4 


1 .5 


0.87 


o 

3. / 


C 1 

5.1 


17 


17 


00 
22 


A C 

45 


00 


AO 

93 


95 


OA 

3u 


IS i jU fJAJ 


e~ 


rliU 


Zz 


oy 
Z5 


A 1 1 

U. 1 i 


A /IQ 
U.4J 


1.0 


c 
J. / 


A 
6.2 


y 


00 
22 


oo 


cc 

JO 


00 
y2 


CO 

j2 






15% 


25 


c o 
J.Z 


0.096 


A on 
0.39 


1 

1.6 


c 
J. 2 


c 
5. / 



8.2 


OA 

2U 


00 
3 / 


5 / 


OA 

yo 


AO 

62 






lU / 


oo 
ZZ 


1.1 


All 

U. 1 1 


A /10 


1.0 


< 
J. / 



0.2 



y 


00 
22 


00 

jy 


CQ 

JO 


01 
y 1 


An 
4/ 






HIO 


17 


34 


0.41 


1.7 


2.4 


7.6 


7.8 


1 1 


26 


47 


69 


97 


48 






15% 


20 


6.5 


0.32 


1.8 


2.3 


7.1 


7.3 


10 


25 


45 


67 


96 


57 






lU / 




1.3 


A 'JO 


1.6 


A 
Z.4 


/.o 




' ' 


00 

2 / 


An 
4/ 


oy 


oo 

y / 


A A 
44 


B = 2(JU ^Cj 




HIO 


17 


27 


A 1 /I 

U.14 


A A 
0.4 


1 o 
l.Z 


3.5 


3.9 


6.6 


20 


Apt 

4U 


Ac 
65 


AO 

y / 


A A 

46 








OA 
ZU 


5.1 


A 1 "J 

0.13 


A OO 

0.3o 


1.1 



3.2 




3.6 


j.y 


1 
18 


00 
38 


AO 

63 


AC 

y5 


cc 
55 






lU/ 


1 Q 
16 


1 


A 1 /I 
U.14 


A A 
U.4 


1 o 
l.Z 


c 


A 
4 


0.0 


OA 
2U 


A 1 
41 


AC 
Oj 


OA 

yo 


A 1 
41 






HIO 


17 


27 


0.2 


0.82 


1.2 


3.7 


4.1 


6.7 


20 


40 


65 


97 


45 






15% 


20 


5.2 


0.15 


0.88 


1.2 


3.5 


3.8 


6 


18 


38 


63 


95 


55 






Tm 
lU/ 


1 o 
1 / 


1 


AID 

u.iy 


A Q 


1 o 
l.Z 


c 

J.O 


A 1 
4.1 


A 
0. / 


OA 
2U 


A 1 
41 


AC 

05 


OA 

yo 


A 1 
41 


D OCA .i/^ 


e~ 


rilU 


1 n 

ly 


o< 
25 


A AA/1 

u.uy4 


A OA 

o.zy 


A QA 


A 

2.4 



2.8 


A A 

4.y 


1 

10 


1A 

34 


CQ 

5is 


AC 

95 


An 
4/ 






15% 


21 


4.7 


0.088 


0.27 


0.83 


2.2 


2.5 


4.4 


15 


32 


56 


93 


56 






T07 


19 


0.95 


0.094 


0.29 


0.85 


2.4 


2.8 


5 


16 


35 


58 


95 


42 




e± 


HIO 


18 


26 


0.14 


0.56 


0.88 


2.6 


3 


5.2 


17 


36 


60 


96 


46 






15% 


21 


4.9 


0.11 


0.6 


0.86 


2.4 


2.7 


4.6 


15 


34 


58 


94 


55 






TO? 


18 


0.99 


0.13 


0.55 


0.89 


2.6 


3 


5.2 


17 


36 


60 


96 


42 


B = 300 


e~ 


HIO 


19 


26 


0.069 


0.23 


0.66 


1.7 


2.2 


4.2 


15 


33 


57 


95 


47 






15% 


21 


4.9 


0.063 


0.22 


0.66 


1.6 


1.9 


3.7 


14 


31 


55 


93 


56 






T07 


19 


0.99 


0.068 


0.23 


0.67 


1.8 


2.2 


4.3 


15 


33 


57 


95 


42 




e± 


HIO 


18 


27 


0.089 


0.37 


0.67 


1.8 


2.3 


4.3 


15 


34 


58 


95 


47 






15% 


21 


5 


0.071 


0.39 


0.66 


1.7 


2 


3.8 


14 


32 


56 


94 


56 






T07 


19 


1 


0.085 


0.36 


0.68 


1.8 


2.3 


4.4 


16 


34 


58 


95 


42 


B = 400 /iG 


e~ 


HIO 


18 


33 


0.057 


0.26 


0.93 


3.1 


4.1 


7.5 


24 


48 


72 


97 


48 






15% 


20 


6.2 


0.05 


0.24 


0.9 


2.8 


3.7 


6.7 


22 


46 


71 


95 


57 






T07 


18 


1.2 


0.056 


0.26 


0.94 


3.1 


4.1 


7.6 


25 


49 


72 


96 


43 




e± 


HIO 


18 


33 


0.56 


2.3 


1.5 


4.1 


4.5 


7.7 


25 


48 


72 


96 


48 






15% 


20 


6.2 


0.41 


2.6 


1.5 


3.8 


4.1 


6.9 


23 


46 


70 


95 


57 






T07 


18 


1.2 


0.52 


2.2 


1.5 


4.1 


4.5 


7.8 


25 


49 


72 


96 


43 


Starbursts: NGC 253 Core 


B = 50 /iG 


e 


HIO 


6 


13 


0.22 


0.84 


3.7 


14 


14 


16 


32 


53 


78 


97 


30 






15% 


6.5 


2.4 


0.19 


0.76 


3.4 


13 


13 


15 


30 


50 


76 


98 


36 






T07 


5.8 


0.48 


0.21 


0.85 


3.8 


14 


14 


17 


32 


53 


79 


97 


27 






HIO 


6 


13 


0.96 


3.9 


4.7 


16 


14 


17 


32 


53 


78 


97 


30 






15% 


6.5 


2.4 


0.72 


4.3 


4.3 


15 


13 


16 


31 


51 


76 


98 


36 






T07 


5.8 


0.48 


0.9 


3.8 


4.7 


16 


15 


17 


33 


53 


79 


97 


27 


B = 100 mG 


e" 


HIO 


9.9 


19 


0.13 


0.54 


2.4 


9.5 


10 


14 


31 


54 


80 


100 


45 






15% 


11 


3.6 


0.12 


0.49 


2.2 


8.7 


9.7 


13 


30 


52 


78 


100 


53 






T07 


9.9 


0.72 


0.13 


0.54 


2.4 


9.6 


11 


15 


32 


55 


80 


100 


41 




e± 


HIO 


10 


19 


0.65 


2.7 


3.1 


11 


11 


15 


31 


54 


79 


100 


46 






15% 


11 


3.5 


0.48 


3 


2.9 


9.9 


10 


14 


30 


52 


77 


100 


54 






T07 


10 


0.71 


0.61 


2.6 


3.1 


11 


11 


15 


32 


55 


80 


100 


42 


B = 150 


e~ 


HIO 


20 


21 


0.086 


0.36 


1.5 


5.6 


6.4 


9.2 


22 


40 


62 


94 


53 






15% 


23 


4 


0.077 


0.33 


1.4 


5.2 


5.8 


8.4 


21 


39 


61 


93 


62 






T07 


21 


0.81 


0.085 


0.36 


1.5 


5.6 


6.4 


9.2 


22 


41 


63 


94 


48 




e± 


HIO 


20 


21 


0.39 


1.6 


1.9 


6.4 


6.8 


9.5 


22 


41 


63 


94 


53 






15% 


23 


4 


0.29 


1.8 


1.8 


5.9 


6.2 


8.7 


21 


39 


61 


93 


62 
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TABLE 4 — Continued 



Model 


Prim, /sb'' £„orra" Scale-i 






Maximum fraction of 7-ray background'^ 








100 keV 


300 keV 


1 MeV 


3 MeV 10 MeV 30 MeV 100 MeV 300 MeV 1 GeV 


10 GeV 


100 GeV 




GeV % 


% 


% 


% % % % % % 


% 


% 



B = 200 ^G e 



B = 250 tiG e 



B = 300 tiG e 



B = 400 AtG e 





T07 


21 


0.81 


0.36 


1.6 


1.9 


6.4 


6.8 


9.5 


22 


41 


63 


94 


48 


e 


HIO 


4.3 


21 


0.073 


0.33 


1.5 


6.5 


8.3 


13 


35 


64 


91 


99 


40 




15% 


4.6 


4 


0.065 


0.3 


1.4 


5.9 


7.6 


12 


33 


61 


89 


99 


48 




T07 


4.1 


0.81 


0.073 


0.34 


1.6 


6.6 


8.4 


14 


36 


65 


91 


99 


36 


e± 


HIO 


4.6 


21 


0.78 


2.7 


2.6 


8 


9 


14 


35 


64 


91 


99 


40 




15% 


4.7 


4 


0.64 


3 


2.5 


7.4 


8.2 


13 


33 


61 


89 


99 


48 




T07 


4.2 


0.8 


0.75 


2.6 


2.6 


8 


9.1 


14 


36 


65 


91 


99 


36 




HIO 


2.9 


23 


0.24 


0.57 


1.6 


5.1 


6.4 


11 


34 


65 


93 


97 


36 




15% 


3.2 


4.4 


0.22 


0.55 


1.5 


4.7 


5.8 


10 


31 


62 


91 


97 


43 




T07 


2.8 


0.9 


0.24 


0.57 


1.6 


5.1 


6.5 


12 


35 


66 


93 


97 


32 


e± 


HIO 


3.3 


45 


0.41 


2.2 


1.9 


5.4 


6.1 


10 


31 


61 


91 


97 


46 




15% 


3.5 


7.8 


0.35 


2.3 


1.8 


5 


5.6 


9.3 


28 


59 


89 


98 


53 




T07 


3 


1.8 


0.42 


2.1 


2 


5.5 


6.3 


11 


32 


63 


92 


97 


40 


e" 


HIO 


17 


18 


0.039 


0.18 


0.73 


2.9 


3.9 


7 


21 


43 


68 


97 


46 




15% 


20 


3.4 


0.034 


0.17 


0.69 


2.6 


3.5 


6.3 


20 


41 


66 


96 


55 




T07 


17 


0.69 


0.039 


0.18 


0.73 


2.9 


3.9 


7 


22 


43 


68 


97 


42 


e± 


HIO 


17 


18 


0.2 


0.85 


0.94 


3.3 


4.1 


7.2 


21 


43 


68 


97 


46 




15% 


20 


3.4 


0.15 


0.94 


0.89 


3 


3.7 


6.4 


20 


41 


66 


96 


55 




T07 


17 


0.69 


0.19 


0.82 


0.94 


3.3 


4.2 


7.2 


22 


44 


68 


97 


42 


e" 


HIO 


20 


13 


0.022 


0.1 


0.36 


1.4 


2.1 


4.4 


16 


35 


59 


94 


48 




15% 


22 


2.5 


0.019 


0.094 


0.35 


1.2 


1.8 


3.9 


15 


33 


57 


93 


58 




T07 


20 


0.51 


0.022 


0.1 


0.36 


1.4 


2.1 


4.5 


16 


35 


59 


94 


43 


e± 


HIO 


20 


13 


0.054 


0.23 


0.4 


1.4 


2.1 


4.5 


16 


35 


59 


94 


48 




15% 


22 


2.5 


0.042 


0.25 


0.39 


1.3 


1.9 


3.9 


15 


33 


57 


93 


58 




T07 


20 


0.51 


0.051 


0.23 


0.4 


1.5 


2.1 


4.5 


16 


35 


59 


94 


43 



Denotes whether the primary electrons and positrons in the model are all electrons {e~) or half electrons and half positrons (e ). 

Starburst fraction used in the model. HIO: our low starburst fraction case, with the starburst fraction in the lHopkins et alj|201Ct) semianalytic models; 15%: our fiducial medium 
starburst fraction case, where 15% of the cosmic star formation rate at all z is in starbursts; T07: our high starburst fraction case, using the values in lThompson. Ouataert. & WaxmanI 
(2007j). 

^ The energy where the ratio of the predicted 7-ray background to the observed 7-ray background reaches its maximum. 

Maximum allowed scaling of the predicted 7-ray background for that model that does not oveiproduce the observed 7-ray background; it is simply the reciprocal of the maximum 
ratio of the predicted 7-ray background to the obsei-ved 7-ray background. 

^ Maximum fraction of the 7-ray background at each energy allowed by the spectral shape of each model. It reaches 100% at £norm- 



